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Abstract 

An equation of state of cold nuclear matter with an arbitrary isotopic composition is 
studied within a relativistic mean-field approach with hadron masses and coupling 
constants depending self-consistently on the scalar mean-field. All hadron masses 
decrease universally with the scalar field growth, whereas meson-nucleon coupling 
constants can vary differently. More specifically we focus on two modifications of the 
KVOR model studied previously. One extension of the model (KVORcut) demon¬ 
strates that the equation of state stiffens if the increase of the scalar-field magnitude 
with the density is bounded from above at some value for baryon densities above the 
saturation nuclear density. This can be realized if the nucleon vector-meson coupling 
constant changes rapidly as a function of the scalar field slightly above the desired 
value. The other version of the model (MKVOR) utilizes a smaller value of the nu¬ 
cleon effective mass at the nuclear saturation density and a saturation of the scalar 
field in the isospin asymmetric matter induced by a strong variation of the nucleon 
isovector-meson coupling constant as function of the scalar field. A possibility of hy- 
peronization of the matter in neutron star interiors is incorporated. Our equations of 
state fulfill majority of known empirical constraints including the pressure-density 
constraint from heavy-ion collisions, direct Urea constraint, gravitational-baryon 
mass constraint for the pulsar J0737-3039B, and the constraint on the maximum 
mass of the neutron stars. 


Contents 

1 Introduction 

2 General formalism 

3 Setup of the RMF models without hyperons 
3.1 Traditional models and the KVOR model 


0 

[13 

Hi 

14 


Preprint submitted to Elsevier 


17 March 2016 





3.2 


w-scaling cut models 


3.3 Scaling functions of KVORcut and MKVOR models 

4 Energy per nucleon, Landau parameters and speed of sound 

5 Comparison of the EoS characteristics with experimental constraints 

5.1 Optical potential 

5.2 Low density behavior of EoS 

5.3 Constraint from particle flow in heavy-ion collisions 

5.4 Direct Urea constraint 

5.5 Gravitational mass versus baryon mass constraint 

5.6 Constraints on the maximum neutron star mass and radii 

6 Strangeness content. Hyperons, 4> mesons, aH scaling. 

7 Conclusion 

A Matching the RME EoS with the crust EoS 
B Derivation of Landau parameters 
References 


0 


24 


27 


27 


31 


34 


37 



M 


57 


1 Introduction 


It is convenient to describe the equation of state (EoS) of the baryon mat¬ 
ter at densities relevant for heavy-ion collisions and neutron stars in terms of 
relativistic mean-held (RMF) models, cf. [T1I21I3] and references therein. The 
original Walecka model [1] included interaction of nucleons with mean helds 
of scalar (a) and vector (ca) mesons. Then an isovector (p) meson held was 
included to ht better the symmetry energy coefficient. The coupling constants 
of the meson helds with nucleons were htted to describe the nuclear satu¬ 
ration, namely, the equilibrium density and the corresponding binding and 
symmetry energies. Next, a a held-dependent potential U{a) = ba^/3 + ca‘^/A 
with two extra parameters b and c was added [5], which allowed to ht val¬ 
ues of the nucleon ehective mass and the incompressibility coefficient. In all 
these models the nucleon mass is considered as a function of the a meson 
held, whereas the masses of the scalar, vector and isovector meson helds are 
treated as constants. Coupling constants are also assumed to be density and 
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field independent quantities. At present, there exists a vast number of general¬ 
izations and modifications of the RMF models. They differ by new terms in an 
effective Lagrangian related to new helds and their interactions P7] . Further, 
RMF models with density dependent nucleon-meson coupling constants were 
developed, cf. Refs. [8II9|I10|I11|I12|I13|I14] . A new density functional DD-MEN 
based to a large extent on microscopic ab initio calculations in nuclear matter 
was presented in [15], which includes nucleons interacting with a, ca, p and 
scalar-isovector 5 meson mean helds with density dependent meson-nucleon 
coupling constants. For a detailed description of the models with density de¬ 
pendent coupling constants see Refs. [TB] . 

On the other hand, various experiments indicate a possibility of modihcations 
of hadron masses and widths in the medium (e.g., see Id)- These changes 
might be related to a partial chiral symmetry restoration in dense and/or 
hot nuclear matter, cf. [MIS] . Motivated by the arguments of a partial chiral 
symmetry restoration at high baryon densities, the Brown-Rho hadron mass 
scaling assumption [20], and the equivalence theorem between different RMF 
schemes. Ref. [2l] demonstrated how one can construct RMF models that 
allow to incorporate simultaneously in-medium modihcations of the baryon 
and meson masses and coupling constants. In this approach the coupling con¬ 
stants and meson masses depend on the a meson held. The masses of the cr, 
w, and p helds, as well as the nucleon mass, are scaled by a universal scal¬ 
ing function d), dependent on the a mean held. A support for the common 
dropping of the N, a , u, p masses comes from the lattice QCD in the strong 
coupling limit [22], where it was found that meson masses are approximately 
proportional to the equilibrium value of the chiral condensate. Remarkably, 
in the case of inhnite matter the ehective hadron mass {ml) and the coupling 
constants {ph) enter all relations only in the combination ml^/g\ that greatly 
simplihes consideration. In order to obtain reasonable EoS, in [ 21 ] also the 
meson-nucleon coupling constants were scaled with the a mean held. Diher- 
ences in the scaling functions for the ehective masses of the u and p helds and 
their coupling constants allow to get a stiher EoS at high densities, a larger 
maximum mass of neutron stars and a sufficiently high critical density for the 
direct Urea (DU) neutrino process, as favored by “the nuclear medium cooling 
scenario” of neutron stars, cf. [221I24II25II2B] . 

Recently, Refs. [ 2328 ] introduced the ring-diagram EOS obtained from re¬ 
alistic low-momentum NN interactions Uow-k- Several microscopic NN po¬ 
tentials (CDBonn, Nijmegen, Argonne V18 and BonnA) were employed. It 
was demonstrated that results are largely improved with an inclusion of a 
Skyrme-type three-body force based on the Brown-Rho scaling, when the 
in-medium meson masses, particularly of cr, oj and p, are slightly decreased 
compared to their in-vacuum values. Best hts were done with the Brown-Rho- 
Ericson (BRE) scaling, cf. [29], given by = 1/(1 -|- Dn/noY^^ with 

D ~ 0.35 ± 0.06, n is the baryon density, no is the nuclear saturation density. 
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and rriN is the nucleon mass in vacuum. Coupling constants were not scaled. 
Reference [30] argued on the basis of the renormalization-group analysis that 
while the pNN coupling constant may fall rapidly with the density, the ujNN 
coupling constant should be scaled at a slower pace, being in favor of the as¬ 
sumption of [21] that the scalings of coupling constants are not universal. In 
the renormalization group approach [30] the effective nucleon mass undergoes 
a drop roughly linear in the density up to a value slightly above the nuclear 
matter saturation density and then stays roughly constant up to the dilaton- 
limit hxed point. Below we will demonstrate a similar behavior within one of 
our selected RMF models. 

Reference [31] formulated most important constraints extracted from the em¬ 
pirical data, which should be fulhlled. It was shown that the models suggested 
in [2T] (namely, MW(nu, z = 0.65) and MW(nu, z = 2.9) labeled then in [ 21 ] 
as KVOR and KVR models, respectively) satisfy majority of the constraints 
known to that time. In the KVR model parameters were adjusted to describe 
the microscopic Akmal-Pandharipande-Ravenhall (APR) EoS of the Urbana- 
Argonne group A18 -|- 5^ -|- UIX* [32] (in relativistic HHJ-parameterization 
of [ 23 ]) at densities below four times the saturation density. A slightly mod- 
ihed parameter set (KVOR) of this RMF model allows for a higher maxi¬ 
mum neutron star mass. The scaled hadron mass-couplings (SHMC) model 
of Refs. [MIISS] is a generalization to hnite temperatures of the KVOR model 
used in PI3T] to describe cold nuclear matter. The SHMC model was suc¬ 
cessfully applied in Ref. [5T]I35II36II37] to description of various characteristics 
of heavy-ion collisions in a broad density-temperature region. 

Some of the previously used constraints were recently tightened and new con¬ 
straints were formulated. At present any useful EoS of the cold hadron matter 
should: 

(z) satisfy experimental information extracted from the description of nu¬ 
clear matter at low densities and not to contradict results of microscopically 
based approaches, 

{ii) satisfy experimental information extracted from the description of global 
characteristics of atomic nuclei, reproducing baryon optical potentials at low 
energies, the binding energy per particle, the symmetry energy coefficient 
and their density derivatives at the appropriate value of the saturation nu¬ 
clear density, values of the incompressibility and the effective nucleon mass, 
(in) not to contradict the direct and elliptic flow constraints [ 23 ] and data 
on the K~^ production in heavy-ion collisions [ 23 ], 

{iv) allow for the heaviest known neutron stars PSR J1614-2230 with the 
mass 1.97 ± 0.04Mq and PSR J0348-I-0432 with the mass 2.01 ± O.O4M0 
[401ITT] (Mq is the solar mass, Mq = 1.99 x 10^^ g), 

(n) avoid a too rapid cooling of the majority of the known pulsars by the 
direct Urea (DU) neutrino processes n—p + e^n-|-z/e [ 23131 ], 
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{vi) explain the gravitational mass and total baryon nnmber of pnlsar PSR 
J0737-3039(B) with at most 1% deviation from the baryon nnmber predicted 
for this particular object [12], 

{vii) yield an appropriate mass-radius relation |13]I44II45| . 

Additionally, the so-constrained EoS of the cold nuclear matter with an arbi¬ 
trary isospin composition should be extended to non-zero temperatures to be 
tested against a wealth of data extracted from analyses of heavy-ion collisions. 

Further information about the compact star mass-radius relation can be found 
in Refs. [ISIHT] . The statistics of 61 measured masses of neutron stars in bi¬ 
nary pulsar systems is collected in [IS]. In contrast to mass determinations, 
there are no high-accuracy radius measurements. In Ref. [13] the neutron star 
radius was constrained to be larger than 11.1 km at the 3a conhdence level, 
with the assumption that the mass of the object is M = 1.76M0, for all com¬ 
binations of other parameters. If taken in combination with the best available 
mass measurement from radio timing [12], the derived constraints prove to 
be inconsistent with all EoSs except very stiff ones. It is interesting to note 
that the lower bounds on the conhdence intervals are still tighter. The X-ray 
spin-phase-resolved spectroscopic study of three thermally emitting isolated 
neutron stars and hts of highly magnetized atmospheric models [H] allowed 
to estimate their compactness, which would also prefer a stiff EoS. None of 
EoSs analyzed in Ref. |3T] is able to explain the suggested radius i? > 14 km 
of the source RX J1856.5-3754. Note that [16] suggested a smaller most prob¬ 
able radius of RX J1856.5-3754, which would be then compatible with some 
known EoSs. A possible correlation between the neutron skin of and the 
neutron star radius was explored in Ref. [50] . 

If we denote as Mp’u the critical mass of the neutron star, at which the nucleon 
DU reactions n -^p + e + 9e, P + ^ n + Ve become possible, every star 
with a mass only slightly heavier than will be efficiently cooled by the 
DU-processes even in the presence of superfluidity (except the case where 
a proton pairing gap is artificially increased) and becomes almost invisible 
for thermal detection within a few years [23E1]. The Il-supernova explosion 
scenario ]5T] and the population synthesis models, see [52] , suggest that most 
of single neutron stars have, probably, masses below 1.5Mq. Therefore one may 
assume that majority of the pulsars, which surface temperatures have been 
measured, have masses M < l.bM©. Thus, [31] suggested as a “strong” DU 
constraint that should be larger than 1.5Mq and as a “weak” constraint 
that at least > 1.35Mq, where 1.35M0 is taken as the mean value 

of the neutron star mass measured in binaries. One can not finally exclude 
a possibility of lower values of (even < I. 35 M 0 ) within a more 

exotic explanation of the present cooling data. However, the absence of the 
DU process for typical neutron star configurations can be considered as the 
most realistic scenario. Reference [2T] demonstrated that the standard RMF- 
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based models yield very low values for the DU threshold density and the 
corresponding star mass, whereas with the help of the cr-dependent scaling of 
hadron masses and coupling constants the problem can be avoided. Similar 
results are obtained in some RMF models with density dependent couplings, 
see [3T] . 


The nuclear EoSs of various models can be characterized by comparing the 
parameters in the expansion of the energy per nucleon near the nuclear sat¬ 
uration density no in terms of the density deviation e = (n — no)/no and the 
asymmetry [3 = (n„ — n^jn for small /3: 


1^) — + + 0{I3^) , 




= J + /^+ 

^=0 O to 


( 1 ) 

( 2 ) 


In this form the EoS is characterized at saturation by the energy per nu¬ 
cleon So, the incompressibility K, the skewness parameter K for the isospin- 
symmetric part and by the symmetry energy J, the symmetry pressure L, and 
the curvature parameter iFsym for the isospin-dependent part of Ssymin). 


The parameters characterizing nuclear matter properties at the nuclear satura¬ 
tion are only partially constrained from the data and have rather broad uncer¬ 
tainties, cf. [53]- One usually hxes the saturation density to be no — 0.16±0.015 
fm“^ and the energy per nucleon. So — —15.6 ± 0.6 MeV. Experimental value 
of the nuclear incompressibility at nuclear saturation is usually extracted from 
the data on the Giant Monopole Resonance (GMR), see [5l]. Reference [55] 
adopted the value K = 240 ± 20MeV, whereas based on the most precise up- 
to-date data on GMR energies of Sn and Gd isotopes, together with a selected 
set of data from ^®Ni to ^°®Pb, Ref. [56] extracted a value of K varied in a 
broad range from 250 to 315 MeV. As has been demonstrated in Ref. m, the 
maximum value of the neutron star mass is a sharp function of the chosen value 
of the (Dirac) effective nucleon mass at the nuclear saturation m^(no) and a 
smoother function of K. The effective nucleon mass is not well constrained 
from the data. For example, in order to describe the spin-orbit splitting in 
atomic nuclei, a value below 0.64m7v is required, see Larger values of the 
effective Dirac nucleon mass are motivated by chiral effective theory expan¬ 
sions |58| and W htting the single nucleon spectra in nuclei [59] with a large 
Landau mass [IJ — 0-9—l.OmAr. The works [60] hnd the Landau mass 
—0.82mjv from the analysis of neutron scattering off lead nuclei. 

^ In our relativistic mean field approximation the Dirac effective mass is the mass 
which enters the Dirac equation for fermions in the medium, whereas the Landau 

mass follows from equality pf/ut-l = (de(p)/dp)|pF, where e(p) = -|- V, 

V is the contribution from a vector potential. We will denote the nucleon Dirac 
effective mass simply as mjy and the nucleon Landau mass as 
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An extrapolation from finite nuclei to nuclear matter of the nucleon optical 
potential at a low particle energy [6T] is better fitted for m*^ ~ 0.7mAr. The 
latter values relate to the Dirac mass m*^ ~ 0.7—O-Sm^v- 

Note that the symmetry energy is often determined as 

^sym(^) 0) Sijlyi ^p) • (3) 

Usually the so-de£ned quantity differs only a little from the value Ssyra{n) 
introduced in Eq. ([2]). 

The symmetry energy is critical for understanding the structure of rare iso¬ 
topes, heavy-ion reactions and for many issues in astrophysics including ques¬ 
tions on the DU threshold, and, as we will demonstrate below, the baryon- 
gravitational mass relation. Experimental data on the density dependence of 
the symmetry energy are available mostly for reactions with stable beams |62] . 
Forthcoming experiments with more neutron-rich nuclei at several radioactive 
beam facilities under construction will be able to improve the situation. Exist¬ 
ing constraints on the values of L and i^sym extracted in analyses of different 
experiments are rather controversial. Reference [53] examined the results of 
laboratory experiments that provided constraints on the nuclear symmetry 
energy for the nucleon density n < no- Some of these constraints have been 
derived from the properties of nuclei, others from the nuclear response to 
hadronic and electroweak probes. The analyses of the isospin diffusion and 
the 'n~ j'n^ ratio in heavy-ion collisions making use of the IBUU04 transport 
model constraining the symmetry energy at sub- and supra-sat nr at ion densi¬ 
ties are discussed in [621IM] . Some constraints on the density dependence of 
the symmetry energy (L = 106 ± 46 MeV and A’sym = 127 ± 290 MeV) were 
extracted from comparing predictions for the neutron-proton elliptic-flow dif¬ 
ference npEFD and npEFR with the FOPI-LAND experimental data |65j . 
Analyzing available experimental information on neutron skins measured for 
26 stable nuclei, from '^^Ca to ^^®U, in anti-protonic atoms. Ref. [66] suggested 
a constraint that L = 55 ± 25 MeV. Studying the measurements of asymme¬ 
try skins in nuclei. Ref. [67] arrived at constraints J = (30.2 — 33.7) MeV and 
L = (35 — 70) MeV. The htting of the data on the atomic masses within the 
Hartree-Fock-Bogoliubov nuclear mass models, based on 16-parameter gen¬ 
eralized Brussels’ Skyrme forces [68] favours a value of ~ 30 MeV for the 
symmetry coefficient J. The neutron skin thickness was extracted from the 
electric dipole response of ^ospp jggj^ which allowed to better constraint the 
L-J correlation suggested in cni Studying broadly different mean held inter¬ 
actions, Ref. EH found an empirical relation for the symmetry energy at the 
saturation density, the slope parameter L, and the curvature parameter iPsym- 
Description of the symmetry energy within quantum statistical approach with 
account for the formation of clusters allowed [72] to join the low-density limit 
with quasiparticle approaches valid near the saturation density. A correlation 


7 



between the symmetry energy and the neutron star radius was studied in [73] 
using quantum Monte Carlo methods. Imprint of the symmetry energy on the 
inner crust and strangeness content of neutron stars was analyzed in [71]. 

At low densities the nuclear EoS is rather well constrained from microscopic 
approaches starting with realistic vacuum NN potentials. One may use the 
variational microscopic (APR) EoS based on the Al8+5u+UIX* NN forces [32] 
or the results of relativistic Dirac-Brueckner-Hartree-Fock calculations based 
on the Bonn potentials, see [S] , Monte Carlo simulations [73] or the calculations 
within a chiral effective held theory. Microscopic neutron matter calculations 
based on chiral NN and 3N interactions were confronted with the constraints 
on properties of nuclear matter in [751176] . 

A set of constraints was also formed from a comparative analysis of 240 non- 
relativistic Skyrme parameterizations m describing nuclear matter at densi¬ 
ties n < 3no. Note that Skyrme functionals predict rather strong deformations 
for the light Pb isotopes. According to the analysis of Ref. [78], this contra¬ 
dicts to experimental data on charge radii and magnetic moments of odd Pb 
isotopes. Reference [73] exploited the phenomenological Fayans functional [73] 
with the volume part of the energy htted to the old Friedman-Pandharipande 
parameterization [80], which is rather close to the parametrization of the APR 
EoS. A large set of 263 RMF models was examined in [81] against several sets 
of constraints. Models with held-dependent hadron masses and coupling con¬ 
stants were not considered. 

From experimental results on the single-particle energies of many A-hypernuclei, 
the ehective binding potential acting on A in inh nite nuclear matter with 
the strength Raat ~ —30MeV was extracted in [82]. The S-hyperon potential 
proves to be Uy,n ~ +30 MeV from the observed quasi-free S production spec¬ 
tra. From analyses of the twin hypernuclear formation [33] and S production 
spectra [31] the potential of S was estimated to be around f/=jv ~ —15 MeV. 
With an increase of the nucleon density towards the neutron star interior the 
nucleon chemical potential may become as large as the ehective hyperon mass. 
If so, new Fermi seas of the hyperons should begin to hll in [3511861187] . As the 
result, the EoS becomes softer that causes a decrease of the maximum mass of 
the neutron star with hyperons. The latter may come into contradiction with 
observations of the most massive neutron stars. This is one part of the ’’hy¬ 
peron puzzle”, see [331139] . The puzzle was extensively investigated within the 
RMF models, see in [31190] . Most of the works predicted masses < I. 6 M 0 for 
non-rotating stars PEZ] . Masses < 1.8Mq were obtained in non-relativistic 
phenomenological models [33]. Fully microscopic models based on hyperon- 
nucleon potentials, which include the repulsive three-body forces, also predict 
maximum masses for the neutron stars below the measured values [3T1135] . A 
repulsive three-body hyperon-nucleon force is needed to reproduce the ground 
state properties of medium-light A hypernuclei [33]. Reference [37] demon- 









strated that the combined repulsive interactions mediated by the oj and 0 
mesons cause the hyperons to appear only at high densities. The chiral quark- 
meson coupling model [HS] exploited that the baryon structure variation in 
matter is reflected in the a-fleld dependence of the g*^ effective coupling pa¬ 
rameterized in the linear form [99] . The resulting decrease of g*^ effective cou¬ 
pling with the density increase allowed to describe neutron stars with masses 
up to 1.95M0. Reference |100] used the RMF model with a density-dependent 
parametrization of the nucleon-meson couplings of Ref. [13] to constrain den¬ 
sity dependence of the hyperon-nucleon coupling constants. To satisfy the 
requirement that the lower bound on the maximum mass of a compact star is 
2Mq Ref. [101] derived constraints on the a A and crS coupling constants. Ref¬ 
erences [1021110311104] used a non-linear model involving two additional hidden 
strangeness scalar a* and vector 0 mesons coupled to hyperons and quar- 
tic terms involving vector meson fields. At this price the authors fulfill the 
constraint on the maximum neutron star mass M^ax > 2Mq. 

Another part of the “hyperon puzzle” is that the critical densities, at which the 
first hyperon species occurs in the neutron star matter, are low in all models 
(~ 3no), if one uses coupling constants constrained by the SU(6) symmetry, 
see |lU5lllU6j . and when the hyperons arise, efficient DU reactions on hyper¬ 
ons, e.g. A ^ p + e + Ue, P + e ^ A + Ue, are switched on. Thus, there may 
appear a problem with fulfillment of the “strong” DU constraint. The prob¬ 
lem might be resolved, if one exploits large hyperon pairing gaps or imposes 
more general SU(3) symmetry relations instead of SU(6) ones [97lllU7j . since 
then with the help of an extra fitting parameter one may decrease values of 
attractive coupling constants. Below we will suggest another solution of the 
problem performed within SU(6) symmetry approach. 

In this paper our main goal is to construct some effective models of the EoS of 
the cold baryon matter that would incorporate the decrease of hadron masses 
and the change of coupling constants with an increase of the baryon density 
and, simultaneously, would fulfill various constraints known from analyses of 
atomic nuclei, heavy-ion collisions and neutron stars. The work is organized 
as follows. In Sect. |2] we formulate a general formalism describing our RMF 
models allowing for scaling of the held mass-terms and coupling constants with 
a held variable / oc uXaNicr), where Xo-Ar(cr) is the scaling of the aN coupling 
constant. In Sect. [3] we formulate several specific models without an inclusion 
of the strange content. More specifically we consider two types of modifications 
of the previously studied KVOR model. In one type of models (KVORcut) the 
EoS stiffens provided the function f{n) stops growing above some chosen den¬ 
sity larger than the saturation density. Other model (MKVOR) demonstrates 
the stiffening of the EoS in the /3-equilibrated neutron star matter while keep¬ 
ing it rather soft in the isospin symmetric matter (ISM). This is achieved, on 
one hand, by the proper choice of the scaling functions for u and a meson 
contributions and, on the other hand, by the quenching of the scalar field 
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growth in the asymmetric nuclear matter induced by the scaling of the p me¬ 
son contribution. Some results of the MKVOR model were briefly reported 
in Ref. |lU8j . In Sect. 0] we describe some equilibrium characteristics of our 
EoSs. In Sect. 0] we demonstrate how our EoSs fulhll various experimental 
constraints. Then in Sect. [6] we account for a possible appearance of hyperons. 
Interaction with c^-meson mean held is included. In Conclusion we formulate 
our results. Matching of the RMF EoSs and the crust EoS is considered in 
Appendix Cumbersome expressions for Landau parameters are deferred to 
Appendix [HI Throughout the paper we use units h = c = 1. 


2 General formalism 


We start with a generalized RMF model, where effective hadron masses and 
coupling constants are assumed to be a held-dependent [21]. The Lagrangian 
density is given by the sum of baryon £bar, meson mean-held £mes, and lepton 
Cl terms: 


— -^b + -^mes + Cl . (4) 

We include the full ground-state octet of baryons b with nucleons N = (p, n) 
(proton p and neutron n) and hyperons H = (A°, with the param¬ 

eters listed in Table H] 


Df 




T '^Qujb'X.ujb^fi T '^dpbXpbtPii T i 9 <j}b\<f)b4^^ 


M • 


(5) 


Here T;, are the baryon bispinor, denote the Dirac y-matrices and t is the 
baryon isospin operator. Greek indices run over 0,1, 2, 3, Latin indices run over 
1, 2, 3. Within the approach of Ref. [ 21 ] the coupling constants g^i,, gujb, Qpb, 9(f>b 
are made dependent of the a held with the help of scaling functions as gabX^b, 
QijbXojb, QpbXpb, gtpbXtjib- Amoug the meson mean helds we include the scalar 
meson and vector mesons, m = {a,u}, where the 0 meson is operative, if 
hyperons are incorporated into the scheme, 

Cmes = \d^ad^a - -U{a)- 

^pv dv^p dp^vi Ppv duPp dpPu 5 ^pv du4^p dp(pu ■ ( 6 ) 


By omitting, for simplicity, the pp interaction term, we disregard a possibility 
of the charged p meson condensation considered in [21]. Here U{a) is a self¬ 
interaction term of the a held. The bare masses of baryons, mi,, and mesons. 
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Table 1 

Particle properties. 



n 

P 

AO 

s- 

SO 

S+ 

[I] 

1 

[I] 

0 

OS 

1 

h 

m[MeV] 

938 

938 

1116 

1193 

1193 

1193 

1318 1318 0.5 

105 

h 

1 

2 

1 

2 

0 

-1 

0 

1 

1 1 

2 2 

- 


rrim, are replaced in Eqs. ([5]) and (E]) by the effective masses ml = 
m*^ = mm^rn- The Lagrangian 

Cl = (V 

i 

contains contributions of the light leptons, electrons and muons, Z = e, p, and 
'ipi stands here for the lepton bispinor. The lepton masses mi are the bare 
masses, see Tabled! 

We will assume that the a field dependence enters the scaling function Xmh 
and ‘hfe(m) through an auxiliary variable 


/ — gaNXcrN^Cr) - 

rriN 


( 8 ) 


Following Ref. [21] we exploit the universal scaling for the field-mass functions 




( 9 ) 


We suppose that 


= XlonU) , XphU) = XpnU) , ( 10 ) 

and express the hyperon-mass scaling function through the nucleon-mass 
scaling function as 

= ^N{g<THX<TH -) = ^NixaniaH - — f) , , (H) 

^ niH'' ^ ruH ' XoN 


where ^^,11 is a function of /. 

Taking into account the equations of motion for vector fields, the energy den¬ 
sity of the cold infinite matter with an arbitrary particle composition is recov- 
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ered from the Lagrangian of the model in the standard way, see Ref. [21]: 


= Zl^kin("l6<hb(/),PF,6) + ^kin("loPF,0 

b l=e,fi 

m%p . a', n 2 


+ 


-vAf) + 




2^2 2m%r]Af) ^ 

E]^iAm,pp) = (^pp \Jm? + p^ {iP + 2pp) - m^arcsinh(pF/'m)^ , (12) 


where the Fermi momentum of species i is pF,i = {SA riiY^^ ^^id we introduce 
the dimensionless coupling constants Cm = Qaint^n/i^m, M = {a,u,p), and 
the ratios of coupling constants xmh = QmhIqmn- Since g^N = 0 we have to 
dehne the ratios x^b through ^f^Ar: x^b = g 4 >B/gu)N, the 0-£eld contribution 
to the energy density is determined by the constant Here 

we take = 783 MeV, = 1020 MeV. Masses of other mesons and their 
coupling constants enter the energy density only via the combinations Cm- 
The scaling functions X(^ and $p, Xp for vector meson helds enter only in 
the combinations 


vAf) 


€(/) 

XlNif) ’ 


Vpif) 


Kif) 

Xlwif) ' 


(13) 


Therefore, we actually do not need to determine Xuj and <Fp, Xp separately, 
but only p{f) combinations. The potential U{a) of the a self-interaction might 
be hidden in the dehnition of Pa{f) as 


vAf) 


€[cr(/)] 

xInAA)] 


+ 



U[a{f )]. 


(14) 


The equation of motion for the remaining field variable / follows from the 
minimization of the energy density = 0, which in the explicit form 

reads 


TTl^ f 

-^vAf) = -Y.^b^bif)ps{mb<^b{f),PF,b) -ruNpAfA^i}) ■ ( 15 ) 

^o- b 

Here the hrst term on the r.h.s is the standard source of the scalar held 
associated with the baryon scalar density ps- 


Ps{m,pp) 


27r2 


(jnpp \J m? + Pp — m^arcsinh(pF/m)^ . 


(16) 
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The second term arises due to the scaling functions rji, 




2^2 


2m^h5(/) 




CpV'pif) ^ ^2 


Solving Eq. flT^ we determine the mean scalar field / as a function of baryon 
densities. 


The coupling constants of hyperons to vector mesons can be related to those 
of nucleons with the help of the SU(6) symmetry relations [105] : 

2 

gujA — QujT, — 25'a;H — -QuN , QpA — 0 , QpY, — ^QpE — ^Qpjq , 

2\/2 

25 ' 0 a = 25 ' 0 s = g<j>'E = — ^gcoN, gtfiN = o. (18) 

One could, of course, use more general relations dictated by the SU(3) sym¬ 
metry nog, which contain one additional parameter, but we prefer a more 
traditional approach here. The extension is though straightforward. 

In our approach the scalar meson coupling constants are constrained by hy- 
peron binding energies in the ISM at saturation, which are deduced from 
extrapolation of hypernucleus data. For the given hyperon binding energy per 
nucleon the aH and uH coupling constants must be correlated in our 

model as 


(J‘2 

^Wnd(^o) = -^x^nno - x„h ^anifo) [mN - m%{no)] , 


(19) 


where /o is the solution of Eq. (1T5|1 in the ISM at saturation, Up = Un = 
We use the results of Refs. 


'^Wnd(’2-o) = -28 MeV, £^bind(^o) = 30 MeV, £^bind(’2-o) = -15 MeV . 


In our model the composition of the neutron star core is determined by condi¬ 
tions of the /^-equilibrium imposing the relation among the particle chemical 
potentials 




f^n Qi l^e 7 t^i 


A 

drii 


E[f,{ni }], 


i = b, 


( 20 ) 
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and by the electro-neutrality condition 

= ( 21 ) 

i 

where Qi is the charge of the particle of type i. Solving Eqs. (120]) and fl^ 
together one can express the particle densities n* through the total baryon 
density n = 

Once all particle densities and the chemical potentials are known the pres¬ 
sure of the matter in the /^-equilibrium (BEM) can be calculated from the 
expression 


P[n] = ^^ini-E[f{n),{ni}]. (22) 

i 


For densities n < 0.7no we match the RMF EoS with the crust EoS, as de¬ 
scribed in Appendix The hnal neutron star conhguration follows from the 
solution of the Tolman-Oppenheimer-Volkoff equation. 


3 Setup of the RMF models without hyperons 


In this section we dehne the scaling functions and input parameters for several 
RMF models. First we demonstrate how our generalized energy density func¬ 
tional (IT^ reproduces results of the traditional Walecka (W) model, non-linear 
Walecka (NEW) models and the KVOR model developed in Ref. [21]. Then 
we propose the method how one can make the EoSs in a certain RMF model 
stiffer, restricting a decrease of the nucleon mass with a density increase. The 
mechanism is a purely phenomenological approach in order to obtain models 
which comply with the constraints known from the experimental data. Finally 
we formulate our novel RMF model, labeled as MKVOR, which properties we 
will discuss in details in the subsequent sections. 


3.1 Traditional models and the KVOR model 


The standard W model [1] with inclusion of p mesons can be recovered from 
Eq. (IT^ setting 


vM) = vM) = Vpif) = 1, U{f) = 0, <I>/v(/) = 1 - /, (23) 

= 1 for the meson helds. The hyperon and 0-meson contributions are 
suppressed. The W model can reproduce the saturation phenomenon in the 
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Table 2 

Characteristics of the W, NLW and KVOR EoS at saturation. 


EoS 

fo 

no 

K 

m^(no) 

J 

L 

K’ 

K 

'^'■sym 

[MeV] 

[fm“^] 

[MeV] 

[mw] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

W 

-15.76 

0.193 

546.59 

0.56 

33.67 

109 

-1880 

74 

NLW 

-16 

0.15 

210 

0.85 

29.51 

78 

632 

-20 

KVOR 

-16 

0.16 

275 

0.805 

32 

71 

422 

-85 

Table 3 

Parameters of the W 

, NLW and KVOR models. 






EoS 

02 

^2 


6- lO^ 

c- 10^ 

z 



W 

266.90 

195.70 

54.710 

0 

0 

~ 



NLW 

183.68 

64.545 

100.00 

17.788 

396.74 

- 



KVOR 

179.56 

87.600 

100.64 

7.7346 

3.4462 

0.65 


ISM with the appropriate choice of Eo and uq, but only at cost of a 

too high 


value of the incompressibility and a too small value of the nucleon effective 
mass, see Table |2l The values of the corresponding coupling constants are 
listed in Tabled cf. |111] . 

The mentioned dehcits of the W model can be cured by the inclusion of two 
extra terms of a a held self-interaction |5] 


riM) = riAf) = ^, <!>«(/) = 1-/. U(f) = , (24) 

= Xmn = 1 for the meson helds. Hyperon and 0-meson contributions 
are suppressed. We will denote this model as the non-linear Walecka (NLW) 
model. Two extra parameters b and c allow to choose appropriate values of 
the incompressibility and the effective nucleon mass at no- Below we use the 
parameter set for model 2 from Ref. [112] . see Table [3l Note that they use 
niN = 939 MeV there. The corresponding saturation properties of this EoS 
listed in Table [2] are quite reasonable. With an account for additional soft-pion 
modes this model appropriately describes heavy-ion collision data at collision 
energies < 2GeV/A, see |113] . 

Reference [2T] constructed the RMF model with scaled hadron masses and 
coupling constants such that it matches the APR EoS (in the relativistic HHJ 
parameterization of Ref. [33]) up to n < 4no. To fulhll the DU constraint 
Ref. [21] introduced the scaling function r]p{f) ^ 1. Thereby, the MW(nu, 
z = 2.9) model was constructed, see Eq. (63) in [21], being labeled as the 
KVR model in Ref. [HI]. The idea behind KVOR modihcation of the KVR 
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model was to demonstrate that introducing the additional scaling function 
T]^ 1 one can increase the maximum value of the neutron star mass not 

changing signihcantly the KVR EoS for n < Auq. In Ref. [2T] the model was 
labeled as MW(nu, z = 0.65), see Eq. (58) there. We will further exploit 
notation KVOR of Ref. [3l]. The scaling functions of the KVOR model are as 
follows: 





The saturation properties for this model are given in Table |2j For the chosen 
parameters of the model, cf. Table |3l the scaling function rjp has a pole at 
/ ~ 0.7, which is harmless at vanishing temperatures, since f{n) remains 
smaller than 0.7. However it was shown in Ref. [ 31 ] that at high temperatures 
rip may become singular. To cure this problem Ref. [ 31 ] suggested to exploit 
the Taylor expansion of the function rjp in terms of 1 — keeping Erst ten 
terms. For T = 0 this expansion fully reproduces the KVOR result. In the 
scalar meson sector we use in Eq. (IT^ 



(26) 


Xo-N : 


which is equivalent to the choice used in Ref. ED- Hyperons and (/)-meson 
contributions are disregarded. 


3.2 u-scaling cut models 

The modern experimental data indicate that the value of the maximum neu¬ 
tron star mass should be sufficiently high. It invites us to search for a simple 
way to stiffen the EoS at large densities without altering it for densities less 
than several no. In Ref. [114] we demonstrated that within an RMF model 
the EoS becomes stiffer for n > n* if a growth of the scalar held as a function 
of the density is quenched and the nucleon mass becomes weakly dependent 
on the density for n > n*. In [114] such a quenching was achieved by adding 
to the scalar potential a function of / rapidly rising at / ~ /*, where f* is 
determined by n*. We demonstrate now that the same effect can be reached 
by the appropriate modihcation of the ri^{f) function. We will show that the 
steeper is the change of ri^{f) at / ~ and the smaller value fi^ is taken, the 
larger becomes the maximum neutron star mass computed within the given 
model. 
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Table 4 

Parameters of the KVOR and three KVORcut models with the a; scaling cut. 


EoS 




5- lO^ 

O 

r—1 


(l(jj 

buj 

Went 

266.90 

195.70 

54.710 

0 

0 

0.56 

-0.5 

53.30 

NLWcut04 

183.68 

64.545 

100.00 

17.788 

396.74 

0.4 

-0.5 

50.00 

NLWcutOS 

183.68 

64.545 

100.00 

17.788 

396.74 

0.3 

-0.5 

50.00 

NLWcut02 

183.68 

64.545 

100.00 

17.788 

396.74 

0.25 

-0.5 

50.00 

KVORcut04 

179.56 

87.600 

100.64 

7.7346 

3.4462 

0.454 

-0.5 

55.76 

KVORcut03 

179.56 

87.600 

100.64 

7.7354 

3.4161 

0.365 

-0.5 

46.78 

KVORcut02 

184.26 

87.594 

100.64 

9.9934 

-76.383 

0.249 

-0.2 

74.55 


3.2.1 Went model 


Let us start with the simplest W model and introduce the scaling function 

vUf) = 1 + y[ 1 + tanh(6^(/ - /^))], (27) 

which jumps within a narrow interval \ f — fi^\ ~ -C 1 from unity to the 
value fjuj — ^ + ciw The corresponding parameters are listed in Table 01 As in 
the case of the W model, we put = XaN = XojN = XpN = 1- At the nuclear 
saturation density the Went model produces the same saturation parameters 
as the original W model. Hyperons and 0-meson contributions are excluded. 
We will denote this modification of the W model as the Went model. 

The scalar held / as a function of the nucleon density in the ISM is shown 
in Fig. [T](left panel) for the W and Went models. As we see from this hgure, 
the f{n) dependence changes sharply in the Went model for n > l.luo and 
/(n) saturates at the value / ~ 0.52. The neutron star masses as functions of 
the neutron star central density for these two models are shown on the right 
panel in Fig. [H The value of the maximum mass for the Went model, = 

2.78Mq, proves to be higher than that for the W model, Mm»x = 2.56 Mq. 
The DU threshold density and mass, = 2.42 and = 2.64Mq, are 
also increased for the Went model compared to the corresponding values, 
nou = 1-78, Mqu = 1.57M0, for the original W model. 


3.2.2 NLWcut model 

Now we apply the same scaling fl271) to the NLW model. We consider three 
modifications NLWcut04, NLWcutOS, and NLWcut02 with the parameters 
listed in Table 01 
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Fig. 1. Left panel: The scalar field / as a function of the nucleon density in the 
ISM for the W model with parameters in Table [3] and for the Went model defined 
by Eq. (1271) . Right panel: Neutron star mass versus the central density for the 
W and Went models. Dashes mark the threshold star mass and density for the 
beginning of the DU reactions. Bold dots indicate the maximum neutron star masses. 
The horizontal band shows the uncertainty range for the mass of PSR J0348+0432 
(2.01 ±0.04M©). 



Fig. 2. Left panel: Scalar field / as a function of the nucleon density in the ISM 
for the NLW and several NLWcut models with parameters given in Table [H Right 
panel: Neutron star masses versus the central density for the NLW and NLWcut 
models. Dashes denote the thresholds of the DU reactions. Notations of the models 
are given in the plot legend. 

The scalar field / as a function of the density is shown in Fig. [2](left) for the 
original NLW model and the NLWcut models. As we see, the f{n) dependence 
changes sharply for the NLWcutOd, NLWcutOS and NLWcut02 models for n > 
3.8no, 2?7 ,o, 1.4no, respectively, and at higher densities f{n) remains constant 
at values 0.36, 0.27, and 0.22. The maximum neutron star masses for these 
models are shown as functions of the central density on the right panel in 
Fig. [21 As in the case of the Went model, the values of the maximum masses 
for the NLWcut models prove to be higher than that for the original NLW 
model (Mmax = 1-73 M©). The maximum mass increases with a decrease of f^. 
The curve M(ncen) for the NLWcut02 model enters the band 2.01 ± 0.04 M©, 
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corresponding to the pulsar PSR J0348+0432. The masses corresponding to 
the DU thresholds also increase with a decrease of but remain too low to 
satisfy the even weaker DU constraint, > 1.35 Mq. 

We discuss the Went and NLWcut models only to demonstrate how with the 
help of the simple cj-cut procedure one is able to increase easily the maximum 
neutron star mass and the threshold DU density. However, the models fail to 
fulhll some other constraints. For example, the Went model does not pass the 
constraint from the particle flow in heavy-ion collisions. All Went and NWLcut 
models do not pass the constraint on the baryon mass vs. the gravitational 
mass of a neutron star. Therefore now we will try to apply the cut procedure 
to construct a more realistic RMF model. 


3.2.3 KVORcut and MKVOR models 

We construct now a set of new KVORcut models which reproduce the EoS of 
the KVOR model up to a certain value of the density and yield stiffer EoSs at 
higher densities. To obtain such a behavior we modify only the scaling function 
of the u meson for / > /^ as follows 


vM) 


f 1 + ^/o \ 

vi+^/y 


+ 


2 L 


1 -h tanh(6^(/ - /^)) 


(28) 


We consider three choices for f^: models KVORcut02, KVORcut03, and KVOR- 
cut04. The parameters of the models are listed in Table 01 Small changes in 
the parameters appear only, if L is taken too close to fn = fino)- We keep 
z = 0.65 for all KVORcut models. 


The input characteristics of the EoS at saturation (£^bind(^o )5 m^(no)/mjv, 
K, and J) are the same as in the original KVOR model. The other charac¬ 
teristics change the stronger, the closer to /(uq) the parameter is chosen. 
So, only one parameter changes slightly in the KVORcut04 model having 
Ksym — —87 MeV and in the KVORcut03 model having iFsym — —88 MeV, 
whereas iFgym — —86 MeV in the KVOR model. In the KVORcut02 model 
two parameters change with iFgym — —78 MeV and K' ~ —869 MeV. 

The KVOR model uses a sufficiently high value of the effective nucleon mass 
at nuclear saturation m*(no) = 0.805mjv- In Ref. [21] it was shown that the 
smaller m*(no) is in an RMF model, the larger the value of the maximum 
neutron star mass becomes. Thus, another way to increase M^ax is to choose a 
smaller value of m*(no). We formulate now the new model, labeled as MKVOR, 
which combines several mechanisms to stiffen the EoS discussed above. In the 
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Table 5 


Parameters of the MKVOR model and characteristics of the EoS at saturation. 



fi2 

^p 

6- 10^ 

c- 10^ 

d 

a 

2 ; 

du 

boj 

234.15 

134.88 

81.842 

4.6750 

-2.9742 

-0.5 

0.4 

0.65 

0.11 

7.1 

f. 

/3 

7 

fp 

(0) 

Up 

a^p^ 

(2) 

Up 

(3) 

Up 

dp 

Cp 

0.9 

3.11 

28.4 

0.522 

0.448 

-0.614 

3 

0.8 

-4 

6 


fo 

no 

K 

m*f^{no) 

J 

L 

K’ 

K 

'■sym 



[MeV] 


[MeV] 

[m-Ar] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 



-16 

0.16 

240 

0.73 

30 

41 

557 

-159 



MKVOR model the scaling functions are taken in Eq. fll 2 l) as follows: 


vM) = 

vM) = ( 

Vpif) = a 


1 - + yf 


Cl = c - -cy 


1 + ^/o 
1 + zf 

+ + 


+ Y [1 + tanh( 6 i^(/ - f^))] , 


( 0 ) 

p 




1 + afp 


+ (3 exp -7 


1 + dpif ~ fo) + ep{f — fo) 


(29) 



As in the previously discussed models the hyperon and 0 meson terms are 
omitted. Parameters of the model are listed in Table O 


The form and the parameters of the scaling functions in MKVOR model are 
tuned to satisfy best the experimental constraints, that we demonstrate below, 
and to keep a connection to the KVOR model parametrization. Indeed, the 
first term in is the same as in the KVOR model, the function rj^j and the 
first three terms in rjp are basically the re-parametrization of the functions of 
the KVOR model. The new terms, the second one in 7 ^^ and the last one in 
7 p, are added to control the growth of the scalar field with an increase of the 
density. Other parameters of rjp are also fine-tuned. 


3.3 Scaling functions of K VO Rent and MKVOR models 


The scalar held variables f{n) for the KVOR, three KVORcut and MKVOR 
models are shown in Fig. |3]for the ISM (left panel) and for BEM (right panel). 
In all models, except MKVOR in the ISM, the function f{n) saturates at high 
densities at the values smaller than 1. The dependence of the scaling functions 
rjm on / is illustrated in Fig. HI We see that r]^{f) changes abruptly at a certain 
value of / for the KVORcut models that allows to stabilize the function f{n) 
at a certain value. For the MKVOR model the same mechanism of the / 
stabilization is implemented not in the 7 a;(/) function but in the 7 p(/), which 
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Fig. 3. The scalar field variable / as a function of the nucleon density following as a 
solution of Eq. (fT^ in the ISM (left panel) and the BEM (right panel) for the KVOR, 
KVORcut and MKVOR models. Vertical bars on the right panel show the maximum 
densities reachable in stars with the maximum masses for the corresponding model. 



f 




f f 


Eig. 4. Scaling functions 'r]„,r]i^,r}p as functions of the scalar field / for KVOR, 
KVORcut, and MKVOR models. Vertical and horizontal bars indicate the maximum 
values of / reachable in stars with the maximum masses for the corresponding model. 

drops fast for / > 0.52. The rise of rip{f) for / < 0.5 serves to tame the 
growth of the nuclear symmetry energy with an increasing density and to 
provide, thereby, a sufficiently high DU threshold. Note that in the KVOR 
model the rip{f) diverges at / ~ 0.7 but such a large value of / is not realized 
in the neutron star within the KVOR model. How one can formally avoid this 
divergence without a change of rip{f) in the physical region of parameters was 
demonstrated in Ref. [3l] and explained above. In Fig. |5]we depict r] scaling 
functions as in Fig. IHbut now computed for / = f{n) for the BEM, as follows 
from the solution of Eq. flT^ . 

The scaling functions of the effective nucleon mass as functions of the 
nucleon density for various models are shown in Fig. [H]for the ISM (left panel) 
and the BEM (right panel). According to our assumption the same function 
controls also the variation of a, u and p meson masses. In all models the 
effective hadron masses decrease hrst with the density and then saturate at 
some constant values, except the KVOR model and MKVOR model in the 
ISM. The limiting values of the effective hadron masses sX hq (here pre- 
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Fig. 5. Scaling functions rifj,ri^,r]p as functions of the density following from Eq. (fT^ 
for the BEM. Vertical bars indicate central densities in the stars with the maximum 
masses for the corresponding model. 




Eig. 6. The common scaling function for nucleon and meson masses as a function 
of the density calculated according Eq. Q with /(n) from Eq. (1151) for ISM (left 
panel) and BEM (right panel) for various models. 

sented up to 8 no) in ISM are 0.78, 0.68, 0.59 for KVORcut02, KVORcutOS 
and KVORcutOd models, respectively, and in BEM they are 0.78, 0.68, 0.59, 
0.39 for KVORcut02, KVORcutOS, KVORcutOd and MKVOR models, respec¬ 
tively. Note that the behavior m*^{n) demonstrated in Fig. |6]is qualitatively 
in a line with that expected within the renormalization group approach in 
Ref. [30], where the effective nucleon mass undergoes an almost linear drop 
for densities up to a value above Hq and then stays approximately constant. 


The scaling functions of the coupling constants Xm are shown in Fig. [7] as 
functions of the density in BEM for our models in comparison with those 
in the DD [16] and DD-F [31] models, where hadron coupling constants are 
assumed to be density dependent, but meson masses stay constant, = 1- 
We see that the density dependencies of our scaling functions for the aN, uN 
and pN coupling constants prove to be similar to those exploited in the DD 
and DD-F models. For n ^ (2-3) no in the KVORcut and MKVOR models 
and for all n in the KVOR model the pNN coupling constant decreases more 
rapidly with the density, than the ujNN coupling constant, that is qualitatively 
similar to the behavior discussed in Ref. [3D]- 
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Fig. 7. The scaling functions for the conpling constants Xm, rn = a,oj, p, as functions 
of the total baryon density in the BEM for the KVOR, KVORcut and MKVOR 
models. Solid lines with circles and triangles depict the density dependence of the 
coupling constants in the DD [16j and DD-F m models, respectively. 



Fig. 8. The scalar and vector potentials, Eq. (IHH]) . acting on a nucleon in ISM as 
functions of nucleon density for KVOR and MKVOR potentials in comparison with 
the DBHF calculations [1151116] for the Bonn A potential. 


It is instructive to compare the scalar, Sn, and vector, Viv, potentials acting 
on a nucleon in ISM 


S„(n) = f{n) , W = ^0) 

with the results of the Dirac-Brueckner-Hartree-Fock (DBHF) calculations [11511116] . 
We see that for the MKVOR model the potentials are comparable with those 
follow from the new microscopic calculations mg, but the former rise with 
a density faster at small densities than the potentials calculated in [115] . 
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Fig. 9. Energy per particle as a function of the nucleon density for ISM (left panel) 
and for PNM (right panel). Short dash-dotted line shows the APR EoS [32] and the 
solid line with circles stands for the AEDMC simulations 

4 Energy per nucleon, Landau parameters and speed of sound 


The energy per particle as a function of the nucleon density is shown in Fig. [9] 
for the ISM (left panel) and for the purely neutron matter (PNM) (right panel) 
for the KVOR, KVORcut and MKVOR models in comparison with the micro¬ 
scopic APR EoS [32] and the auxiliary field diffusive Monte Carlo (AFDMC) 
calculations [117] . All considered EoSs demonstrate a similar behavior for 
n < 1.5no and begin to deviate substantially from each other demonstrating 
different stiffness with a density increase. 

The energy of a nucleon in the nuclear matter with proton and neutron particle 
distributions fp{p) and fn{p) follows from the Lagrangian in the mean-field 
approximation ([5]) as 


^i(p) = + (p - U(/))2 + V'ci,i(/), i^n,p, 


(31) 


where vector potentials and currents are defined by the relations 

K),p(/) + Vo,nif) _ f \ + np) 

2 ° m%Uf) ’ 

VoAf) - VoAf) _ _ C^pjnp-nn) 

2 ° ^mlvpif) ’ 

MD + Vnif) ClGn+Jp) Vp{f)-Vn{f) C% - t) 

2 rn%p^{f) ’ 2 Am%pp{f) 

P- Viif) 


f 2 d^p „ ^ [2 d^p 

Hi = I T:-^/i(p), Ji = 


(2vr)' 


P’"'” Ap- Any +<«>?,(/) 


= Mp)- 


(32) 
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The scalar field parameter / follows from Eq. (1151) 


2^2 


vM) 


2 vUf) 
L/ vAf) 


+ + VA’{f)(nn - Up) 


Mf) 


ip)i 




vAf) 

2d^p fp{p) + fnjp) 

\/(p - Vi(f))^ + m%iS>%{f) 


Vpif) 

(33) 


We consider the nucleon distributions with a slightly distorted, non-spherical 
Fermi-surface, so that the nucleon current is non-zero, ji ^ 0. For example we 
can take fi{p) = O^pppl -|- e/3 -|- ap/pp^t ] — p), where e <C 1, |a| -C 1 are 
small parameters. The variation of the nucleon energy with respect to nucleon 
distribution defines the Landau parameters. For momenta close to the Fermi 
surface we can write, cf. [m]. 


Fij{n,n') 


Seiipp^ifi) 
npfjippjfi') ■ 


(34) 


We take here into account that the variation of the density under the variation 
of the particle distribution is Srij = njSfj ~ n^e. 


The Landau parameters can be expressed as matrix elements of the effec¬ 
tive nucleon-nucleon interaction matrix in the particle-hole channel, 

Fij{n’,n) = A){xl)AXj)s , (35) 

where Xi is the isospin spinor of the nucleon of type i, (Xp)a = ^la and {xn)a = 
620 , and the interaction matrix can be written in terms of the Pauli matrices 
acting in the nucleon isospin state, 

n) = F{e) + F'{e) ■ (36) 

Here F and F' are functions of the angle between the directions of the Fermi 
momenta of incoming and outgoing nucleons, cos 6 = {n'n). From Eqs. (I35II36I) 
we find relations, cf. ma 

Fnn = Fpp = F 22,22 = F F F' ^ F^p = F 12,21 = F — F' . (37) 

The Landau parameters are usually expanded in the Legendre polynomials 


Fi = {21 + 1)1 ^F{cose)F{e), f; = {21 + i)J ^F{cose)F\e) . (38) 

The Landau parameters Fq and Fi are important characteristics of the nuclear 
Fermi liquid. For the Walecka model they were first calculated in Ref. m- 
Following |111] . we define the dimensionless Landau parameters /o,i as follows: 

fi = NpFi, i = 0,1, Np = 7 PfFf/ {2A), Pp = {GAn/pF^. (39) 
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Fig. 10. The density dependence of the dimensionless Landau parameters /o and /i 
in the ISM and the PNM, and of the parameters /q and /{ in the ISM. Solid curves 
with filled dots show calculations within the Skyrme BSk21 model |119j . Notations 
of other curves are the same as in previous Figs. 

Np is the density of states at the Fermi surface, Ep = \Jpp + is the 
relativistic Fermi energy of the nucleon. For ISM 7 = 4, for the PNM 7 = 2. 
We will consider here these two cases, since for the nuclear matter of arbitrary 
isotopic composition expressions are more cumbersome. 


The incompressibility K, the square of the sound velocity nf and the chemical 
potential in the ISM are expressed through /o(n) as 


K = 


Ep 


(1 + /o), 




pI 

3Epfi 


(1 + /o), h — 


C^,n 




+ ^Pl + rn*^- (40) 


Note that in the relativistic theory the particle energy at the Fermi surface 
plays the same role as the effective mass in the non-relativistic Fermi-liquid 
theory. It can be expressed in terms of the Landau parameter fp. 


Ep = ^pI + = /i(l + ^/i). (41) 

The symmetry energy is connected with /g in the ISM by the following relation: 

4« = ^(l + /i)- (42) 


For our generalized RMF models with scaled hadron masses and coupling 
constants the Fermi liquid parameters in ISM and in PNM are presented in 
Appendix [Bl 


In Fig. [10] we demonstrate the density dependence of the Landau parameters 
/o and fi computed within the KVOR, KVORcut and MKVOR models. The 
results are shown for the ISM (left panel) and the PNM (right panel) in 
comparison with calculations performed in the framework of the Skyrme model 
BSk21 [119] . We see that the Landau parameters in all models demonstrate 
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Fig. 11. The squared speed of sound as a function of the nucleon density for various 
RMF models in the ISM (left panel) and in the BEM (right panel). 

the qualitatively similar behavior. In the ISM for 5 x 10“^ Hq < n < 0.6 Hq 
we get the scalar Landau parameter /o < —1. In this density interval there 
appears the Pomeranchuk instability with respect to the growth of the long¬ 
wave density fluctuations. In the PNM /o > — 1 and such an instability does 
not arise. 

In Fig. [H] we show the behavior of the squared speed of sound as a function of 
the nucleon density in the ISM (left panel) and in the BEM (right panel). The 
sound velocity demonstrates a monotonous rise with a density increase but 
never reaches the velocity of the light (c = 1). We see also that becomes 
negative in the density interval, where the scalar Landau parameter /o < — 1, 
i.e. in the region of the Pomeranchuk instability, see |12n] . 

In the BEM the squared speed of sound bends slightly (not seen visually on 
the plot) at n ~ O.Sno, which is the point of the muon appearance. Such a 
behavior is typical for third-order phase transitions. 


5 Comparison of the EoS characteristics with experimental con¬ 
straints 


5.1 Optical potential 


The baryon optical potential in the ISM is determined as [sum] 

^opt(^) = - s, {2 m, + Sk), (43) 
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Fig. 12. Left panel: Energy dependence of the nucleon optical potential in the ISM 
at n = no calculated for the KVOR and MKVOR models. Hatched area shows 
the extrapolation from finite nuclei to the nuclear matter [6T]. Right panel; Optical 
potentials of S, A and H for the same models. 

where e is the baryon energy, the scalar and vector potentials acting on the 
given baryon in ISM are defined as 

/^2 

Sb = /) - nib, 14 = x^b 2 • ( 4 ^) 

The expression for the optical potential of antibaryons is obtained by replacing 

H ^ -H. 

The dependence of the nncleon optical potential on the nncleon energy in the 
ISM at n = no is shown in Fig. [12] (left panel). The hatched band is the 
optical potential extracted from the atomic nuclens data |122] and recalcn- 
lated to the case of the inhnite nnclear ISM in Ref. [61]. We see that the 
MKVOR model describes the nucleon optical potential rather appropriately 
for energies £ — ntj^i < 400 MeV. To match the data for higher particle energies 
a momentum dependence of the NN interaction would be required, which 
is not present in the mean-held approach. The KVOR model describes the 
nucleon optical potential for low and high energies better than the MKVOR 
model but does not describe it for intermediate energies. Calculations for the 
KVORcut models do not differ from that for KVOR since for n = hq differ¬ 
ences in the relevant parameters of these models are minor. The isovector part 
of the optical potential — U^ptis) is less constrained by the data and 

we do not discuss it, therefore. 

The energy dependence of the hyperon optical potentials is not yet constrained 
by data. Predictions of our models for hyperons are shown in Fig. [T^fright 
panel). Predictions for the optical potentials Uj^pti^) of antinucleons and an- 
tihyperons in the ISM at n = no are depicted in Fig. [131 Predictions for an¬ 
tiprotons are very important in view of the future experiments at FAIR. For 
the KVOR model the antinucleon optical potential was calculated in Ref. [35] . 












Fig. 13. Energy dependence of the antinucleon optical potential (left panel) and of 
the antihyperon one (right panel) for the ISM at n = no. 

A phenomenological value of the antiproton optical potential is limited within 
the range —(100—350) MeV |123j . The available experimental data from p 
atoms |124j and p scattering off nuclei |125j suggest that the depth of the real 
part of the p-nucleus potential in the interior of a nucleus lies in the range 
-(100—300) MeV. 

5.2 Low density behavior of EoS 


The behavior of the nucleon EoS for n < uq was extensively studied within 
various phenomenological and microscopic models including the chiral effective 
held theory. The EoS derived from our RMF model can be, of course, tuned 
to one of realistic parameterizations of the EoS for n < no and vanishing 
temperature, but this is not our goal in the present work. We just continue to 
exploit simple analytical parameterizations of the scaling factors and 

Pp presented above. 

The energy per particle and the pressure for the PNM are shown in Fig. [14] as 
functions of the nucleon density for n < hq. The curves for all three KVORcut 
models coincide with the KVOR model in this density range. The MKVOR 
curves he within the uncertainty region estimated within the chiral effective 
held theory |^. The pressure of the KVOR model goes beyond the estimated 
uncertainty region. The diherence between the energies per particle for the 
AFDMC EoS |117j and our MKVOR parameterization does not exceed 2 MeV 
for n < ng. In comparison with APR EoS Al8+(5u+UIX* the diherence is a 
bit larger. 

Constraints from Ref. [65] on the density dependence of the symmetry energy 
^^sym(^), cf. Eq. ([2|), were deduced from the comparison of theoretical predic¬ 
tions for the diherence in proton and neutron hows with the experimental 
data of FOPI-LAND collaboration for Au-|-Au collisions at 400 MeV/nucleon. 


29 















n/no n/no 


Fig. 14. Energy per particle (left) and pressure (right) for PNM as functions of the 
nucleon density for n < hq. Overlapping bands are evaluated from chiral effective 
held theory for different NN potentials and include uncertainty estimates due to 
the many-body effects m- AFDMC calculations [117j are shown by empty circles, 
and short dash-dotted line stands for the APR EoS. Solid line presents calculation 
with the KVOR EoS, long dash-dotted line is for MKVOR EoS. 


The constraints are shown in Fig. [15] together with the constraints from the 
chiral effective held theory [76] and from the stndy of the isobaric analog 
states (IAS) [67]. In Fig. [15] we plot the symmetry energy derived for onr 
KVOR, KVORcnt02 and MKVOR models. The cnrves for KVORcntOd and 
KVORcntOS models coincide with the KVOR cnrve. Note that the KVOR and 
KVORcnt models yield J = 32 MeV, whereas the MKVOR model corresponds 
to J = 30 MeV. Becanse of this the MKVOR cnrve lies below the KVOR one 
at uq. The valne £’sym(^) is given by Eq. ([3]). For n < uq the difference be¬ 
tween Ssyra{n) and S{n) is minor. It increases however with a density increase 
for n > Uq and Ssymi^) < Ssymin). For n > Uq cnrve £^sym(’^) for the MKVOR 
model (thin dash-dotted line) deviates slightly from the lower bonndary of 
region estimated from the FOPI-LAND data, whereas the deviation is larger 
for the qnantity Ssym{n) (thick dash-dotted line). The lines for the KVOR and 
KVORcnt models he inside the region. 


The nncleon potential, Uisfin) = Sn + Vn with the scalar S and vector po¬ 
tential V dehned in Eq. (jHj), can also be constrained in the ISM at low 
nncleon densities from measnrements of the kaon prodnction in heavy-ion col¬ 
lisions performed by the Kaon Spectrometer (KaoS) [126] . It is shown on the 
right panel in Fig. [15] together with the cnrves prodnced by onr models, see 
also |M|. We hnd that KVORcnt02 model is too stiff to fnlhll the constraint. 
The KVOR, KVORcntOd and KVORcnt03 models fnlhll the constraint except 
near the border 2no. The model MKVOR fnlhlls the constraint fnlly. 
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Fig. 15. Left panel: The symmetry energy coefficients Esyrain), cf. Eq. ([2]), and 
^symin) given by Eq. ([3]) as functions of the nucleon density calculated for various 
models. The bands are constraints from chiral effective field theory (xEFT) |76j . 
from study of IAS and from the FOPI-LAND experimental data for Au+Au 
collisions at 400 MeV/nucleon [65] . Right panel: The nucleon potential Uiy as a 
function of the nucleon density in the ISM for various models. The curve labeled 
as KaoS shows the upper boundary for the potential deduced |126] from the kaon 
production measured by KaoS collaboration. 



Fig. 16. Pressure as a function of the nucleon density for ISM (left panel) and for 
the PNM (right panel). Double-hatched area is the constraint from the particle flow 
in heavy-ion collisions [38], hatched area is the kaon flow constraint extracted in 
Ref |39j from the analysis of Ref. |127j . Line with bold dots shows the extrapolation 
of the pressure consistent with the GMR, cf. Ref. [39]. On the right panel two 
double-hatched areas show the pressure consistent with the flow data after inclusion 
of the isospin asymmetry terms with stiff and soft density dependencies. 

5.3 Constraint from particle flow in heavy-ion collisions 


The analysis [38] of experimental data on transverse and elliptical flows of 
particles in heavy-ion collisions reveals a correlation among the flow magnitnde 
and the stiffness of the EoS in the ISM. The pressnre as a fnnction of the 
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nucleon density in the ISM is shown in Fig. [T6](left panel). Double-hatched 
band ranging from 2no to about 4.5no corresponds to the so-called particle- 
flow constraint of Ref. |3B]- It rules out very stiff EOSs and puts in a challenge 
the requirement that EoS should be stiff enough to fulhll the constraint on 
the maximum neutron star mass extracted from the measurement of pulsar 
masses in Ref. [401I4T] . In Fig. [16] we show also a restriction on the stiffness of 
the EoS from n ~ 1.2?7,o to n ~ 2 .2nQ extracted in |39] from the analysis of 
the kaon flow in Ref. |127] and the extrapolation of the pressures consistent 
with the GMR data analysis [32] • The KVOR model satishes the requirements 
from heavy-ion collisions. However, the maximum neutron star mass with the 
KVOR EoS is equal to 2.01Mq and fits only marginally the experimental value 
of the mass 2.01±0.04Mq of the pulsar PSR J0348-I-0432 [3T|. Since hyperons 
are not included in the KVOR model and their inclusion may only result 
in a decrease of M^ax, the KVOR model needs a revision. The EoS for the 
KVORcut04 model proves to be only slightly stiffer than that of the KVOR 
one: the pressure curve for this model passes very close to the curve for KVOR 
model. The curve for the KVORcut03 model goes over the upper boundaries 
of the hatched regions. The curve for the MKVOR model passes through the 
hatched regions for n < 4no but escapes it at higher densities. Contrarily, the 
curve for the KVORcut02 model lies above the region for n < d.SriQ but enters 
the double-hatched region for higher n. Thus, we see that the KVORcut03 
and MKVOR models are the most promising models for the simultaneous 
fulhllment of the particle-flow and maximum neutron star mass constraints. 

For the PNM the upper and lower double-hatched regions in Fig. [T6| (right 
panel) correspond to the pressure in the ISM consistent with the flow data 
after inclusion of the pressure from the symmetry energy term, £syra{n), with 
a strong or weak density dependence. All presented curves except the KVOR- 
cut02 one fit the double-hatched region with the strong density dependence of 
£syra{n). The curve for KVORcut02 model lies slightly above the upper bound¬ 
ary of the double-hatched region in the region n < 3no that might indicate 
too high stiffness of this EoS at these densities. 


5.4 Direct Urea constraint 


The DU processes on neutrons, u—)-p-|-e-fr'e,p + e—can occur if 
the Fermi momenta of particles satisfy the inequality pF,n < Pf,p + PF,e- For 
the hyperon free matter this inequality can be rewritten with the help of the 
electro-neutrality condition fl2T]) as 


Hp i 

— > a:Du —- rpi — > 

n i + (i + 2:y^)3 


n. 


X. = 


Ue + n. 


(45) 
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where xdu is the DU-threshold proton fraction. The values of the DU threshold 
density, and the corresponding mass of the star, strongly depend 

on the density dependence of the symmetry energy. 

In Fig. [T7](left panel) we show the proton fraction as a function of the den¬ 
sity for EoSs under consideration for BEM, together with DU threshold ratios 
xdu(^)- We see that the ratio xdu(’^) is only slightly model-dependent. Once 
the DU threshold density is reached in the center of a neutron star, the very 
efficient process of neutrino cooling of a neutron star kicks in. As a result every 
star with a mass only slightly exceeding will be rapidly cooled down by 
the DU processes, n —)■ p + e + Ue, P + e —)■ n + u^, even in the presence of super¬ 
fluidity, and becomes almost invisible for the thermal detection within a few 
years The type Il-supernova explosion scenario |ST] and population 

synthesis models [52] predict that most of single neutron stars have, probably, 
masses below 1.5M0. Therefore, one may assume that the majority of pulsars, 
which surface temperatures are measured, have masses M < I.SMq. The av¬ 
erage mass of neutron stars observed in binaries is 1.35Mq. That motivated 
Ref. [3T| to treat the inequality > I-SMq as a “strong” DU constraint 
and the inequality > 1.35Mq as a “weak” DU constraint on the nu¬ 

clear EoS. The appropriate description [25)126] of the data on the cooling of 
the pulsar in Cassiopea A also requires absence of the DU reactions in the 
star interiors. (Except for the cases where one artificially increases the proton 
gap, which we consider as unrealistic.) The analysis of these data in the ex¬ 
isting cooling scenarios supports the constraint that should be > 1.5Mq, 
see [12811^ . Although one can not still exclude the possibility of lower values 
of Mp’u within a more exotic explanation of the present cooling data, e.g. with 
artificially enhanced proton gaps, the absence of the DU process for typical 
neutron star configurations might be considered as the most probable scenario. 

The DU thresholds and the corresponding star masses for various models are 
collected in Table (6] We see that for all models, which we consider here, the 
values of prove to be fairly above the value 1.5Mq required by the DU 
constraint, although the threshold densities are not so large. The reason 
for this is that the EoSs for our models are sufficiently stiff. To show this in 
Table [6] we also indicate values of the maximum neutron star masses for the 
models under consideration. 

From Fig. [T7] (left panel) we also see that in the density interval uq < n < dug 
the proton fraction of the MKVOR model is smaller than those for the KVOR 
and KVORcut models. For n > dng the proton fraction for MKVOR model 
exceeds that for KVOR but remains smaller that those for the KVORcut 
models. Also we see that the smaller the value of in KVORcut models is 
chosen, the higher the proton fraction is. For n < 0.7no we use the BPS EoS 
and the curves for particle concentrations presented in Fig. [T7] (left panel) 
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Fig. 17. Left panel; Proton fractions in the BEM for various models (thick lines) 
and the thresholds of the DU reaction xdu (corresponding thin lines) as functions 
of the nucleon density. Vertical line indicates the upper border of the crust. Right 
panel: Gravitational mass of a neutron star, Mq, versus the baryon mass, Mb, 
for various models. The double-hatched rectangle is the constraint for the pulsar 
J0737-3039(B) [42]. The two empty rectangles show the variation of the constraint, 
when the assumed loss of the baryon mass during the progenitor-star collapse 
amounts to 0.3% Mq and 1%Mq. Line for the KVORcut04 model coincides with 
that for the KVOR model. 

Table 6 

Maximum neutron star masses and DU thresholds for the KVOR, KVORcut and 
MKVOR modeK_ 


EoS 

[Mq] 

[’^o] 

[Mq] 

KVOR 

2.01 

3.96 

1.77 

KVORcut04 

2.09 

3.38 

1.67 

KVORcut03 

2.17 

2.85 

1.68 

KVORcut02 

2.26 

2.43 

1.84 

MKVOR 

2.33 

3.84 

2.14 


should be replaced by that in the model for the crust. 


5.5 Gravitational mass versus baryon mass constraint 


Reference |12] studying pulsar B in the double pulsar system J0737-3039 sug¬ 
gested a test on the EoS of nuclear matter. The system J0737-3039 consists of 
a 22.7 ms pulsar J0737-3039A, and a 2.77 ms pulsar companion J0737-3039B, 
orbiting near the common center of mass in a slightly eccentric orbit of 2.4 h 
duration. The gravitational mass of the pulsar B is carefully measured in |129] 
to be 1.249 TO.OOIMq. Such a low mass could be an indication that the pulsar 
B was formed in a type-I supernova of an 0-Ne-Mg white dwarf driven hydro- 
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statically unstable by electron captures onto Mg and Ne. Such an instability 
occurs when the progenitor star core density reaches a well-dehned critical 
value (~ 4.5 x 10® g/cm^) which corresponds to a well-dehned mass of ONeMg 
core (~ 1.37 Mq). Assuming that the loss of matter during the formation of 
the neutron star is negligible, Ref. |12] predicted the baryon mass Mb = u Nb 
to be Mb(J0737-3039B) = 1.366 — 1.375Mq, where Nb is the total number 
of baryons in the star and u stands for the atomic mass unit equal to 931.5 
MeV. The baryon number is calculated as 


R 

NB = 4:Tr J 

0 


drr'^nir) 


^1 - 2GAf(r)/r’ 


(46) 


where G is the gravitational constant, R is the radius of the star, n{r) is the 
baryon number density and M(r) is the total mass accumulated inside the 
radius r. 


If the formation mechanism of the PSR J0737-3039 system and the assump¬ 
tion of a negligible baryon loss of companion B during its creation are valid, 
the result of Ref. |12] provides a strong constraint on the nuclear EoS. Micro¬ 
scopically motivated EoSs like the relativistic Dirac-Brueckner-Hartree-Fock 
EoS [nn], the APR EoS Al8+6v+mX* [52] and the AFDMC one [II7|, as 
well as EoSs of many RMF-based models do not fulhll this constraint. The 
baryon loss and variations of the critical mass due to carbon flashes during the 
collapse may result in a lowering of Mb by ^ 1% Mq, that partially softens the 
constraint. However, many EoSs do not satisfy even this weaker constraint; 
see Ref. [5T] . 

In Fig.[T7](right panel) we show the gravitational neutron star mass Mq versus 
its baryon mass Mb- The two empty rectangles show the change of the con¬ 
straint, when the assumed loss of the baryon number in the collapse amounts 
to 0.3%Mq and to 1%Mq. Approximately the same (from O.3%M0 to 1%Mq) 
constraint box was proposed in Ref. |131j . which found in their model that 
the mass loss of the collapsing 0-Ne-Mg core during the explosion leaves the 
neutron star with a baryon mass of M = 1.36 ± O.OO2M0. We see that the 
KVORcut02 EoS does not match even a weak constraint, when the assumed 
baryon loss in the course of the star collapse amounts to 1%Mq. The KVOR- 
cut04 line coincides with that of the KVOR since both EoSs differ only a 
little for n < ricenlMc = I. 25 M 0 ] under consideration. The KVORcut03, 
KVOR and KVORcut04 curves match marginally the weak constraint. The 
MKVOR model hts marginally the “strong” constraint (the curve touches the 
left boundary of the hatched box). Comparing the hgures from the left and 
right panels we can End a correlation: the smaller the proton fraction is within 
the density interval no < n ^ 2.5no, the better the given EoS satisfies the Mq- 
Mb constraint. The KVORcut02 model yields the largest proton fraction in 
the given mass interval and does not fulhll the constraint. The proton fractions 
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Table 7 

Parameters of the scaling function (H71) . the corresponding values of L and ilTsym 
and the central densities (ucen) of the neutron star with the mass 1.2b Mq. 


/3 

7 

L [MeV] 

K,yra [MeV] 

“^cen [^ 0 ] 

0 

0 

85.30 

18.77 

2.33 

0.80 

3.60 

70.72 

-80.90 

2.61 

0.80 

7.50 

54.91 

-155.4 

2.70 

1.50 

7.50 

30.94 

-189.3 

2.68 


for KVOR, KVORcutOd and KVORcutOS models are almost indistinguishable 
for n < 2 no and their Mq-Mb lines almost coincide in the considered inter¬ 
val of masses. The MKVOR model produces the smallest proton fraction for 
n < Ario and it matches marginally even the strong constraint. 

Let us elaborate more on the connection between the proton concentration 
and the Mq-Mb constraint. The proton concentration is determined by the 
symmetry energy and therefore is correlated with the values of J and L in 
Eq. ([T]). The value of J may vary only a little from 28 MeV until 34 MeV 
or even in a more narrow interval, whereas the value of L is less known at 
present [391163] . From the results shown in Fig. [T7]we see that all KVORcut 
models, which miss the Mq-Mb constraint, have a larger value of L than the 
MKVOR model, which satishes the constraint better. To study how much the 
fulhllment of the Mq-Mb constraint is sensitive to changes of L, we consider 
a set of test models with all parameters being the same as for the MKVOR 
model, except the scaling function ? 7 p(/) replaced by the following one 



Varying the parameters and 7 we can simulate different values of L, see 
Table [71 In the upper row in Fig. UHl we plot the corresponding symmetry 
energy and the proton concentration as functions of the density. For a hxed 
density n, the decrease of L in our cases leads to a decrease of Ssym if n > no 
and to an increase of Ssym if n < uq. Correspondingly, the proton concentration 
at a hxed density n decreases with the decrease of L if n > no and increases 
for n < no. The EoS softens for n > no herewith and the central density 
needed to support the neutron star with the hxed mass 1.25 Mq increases 
from 2.33 no to 2.70 no for L decreasing from 85 MeV to 55 MeV. Further, 
for L decreasing from 55 MeV to 24 MeV it decreases weakly from 2.70 no to 
2.66 no. The proton concentration decreases with an L decrease for all r, since 
the radius of the star decreases with decrease of L. 

In Fig. dniwe depict the accumulated gravitational and baryon masses of the 
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Fig. 18. Upper row: The symmetry energy (left panel) and the proton fraction (right 
panel) as functions of the nucleon density. Lower row: The nucleon density and the 
proton fraction as functions of the radial coordinate, r, for the neutron star with the 
mass 1.25 Mq. Calculations are performed in the RMF model with the parameters 
and the scaling functions as in the MKVOR model and the rjp function given in 
Eq. (I47ji with parameters from Table [3 for various values of L. The rounded values 
of L label curves. 

neutron star as functions of the radial coordinate. We observe that the baryon 
mass is larger than the gravitational mass, Msir) > Mdr), for any r. The 
radii at which the gravitational and baryon masses saturate, are almost the 
same and become smaller when the value of L decreases. The right panel in 
Fig. [19] demonstrates that the more we decrease the L value the closer the 
baryon mass approaches the empirical “strong” constraint from [l2] shown by 
the double-hatched band, and enters it for L < 30.9 MeV. 


5.6 Constraints on the maximum neutron star mass and radii 


Recent precise determinations of the masses of pulsar PSR J1614-2230 in 
Ref. [10] as 1.97 ± O.O4M0 and of pulsar PSR J0348-I-0432 in Ref. [H] as 
2.01±0.04 Mq delimit the EoS of neutron star matter. Other masses of the very 
massive pulsars are determined with large experimental error bars and with 
additional theoretical uncertainties. Reference [132] performed initial spectro¬ 
scopic observations and additional photometry from the y-ray pulsar PSR 
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Fig. 19. Accumulated values of gravitational and baryon masses as functions of the 
radial coordinate, r, for various L. Calculations are performed for the same models 
as in Fig. [THl The double-hatched band on the right panel shows the empirical 
constraint from Ref. |42j . 


J1311-3430 with the period 94 min. Simple heated light-curve hts give the 
estimate for the pulsar mass 2.7Mq. Incorporating systematic light-curve dif¬ 
ferences the authors estimated M > 2.1Mq. For further information about 
the compact star mass-radius relation see Refs. [461177] . 

In Fig. [20] we demonstrate the neutron star masses and radii as functions of 
the central density for our EoSs with the crust computed with the very same 
EoS (we name it “without crust”) and with the crust computed with BPS 
EoS (named ’’with crust”). Bold dots show maximum masses. We see that 
the calculation “with crust” affects only a little the values of the neutron star 
masses increasing them not more than by 0.05 Mq and increasing the radii by 
about 1 km for the masses corresponding to central densities 2-4 uq and by 
~ 0.5 km for heavy stars with central densities 5-6 uq. The band shows the 
uncertainty range of the masses for PSR J0348-I-0432 (2.01 ± 0.04Mq). All 
the models fulhll the M^ax constraint. We see that for the KVORcut models 
the lesser the value f^j is chosen, the larger the value of the maximum mass 
is and the smaller the central density corresponding to Mm ax becomes. The 
MKVOR EoS is the stiffest among EoSs considered here. The central density 
Urpu ( Mmax ) for the MKVOR EoS is above that for KVORcut02 EoS. Dashes 
on the right panel indicate the DU thresholds. The MKVOR model yields the 
highest DU threshold mass Mq^, cf. Table El 

The internal structure of the neutron star with the maximum possible mass 
is illustrated in Fig. [21] for the various EoSs we study. In the left column we 
show the adiabatic index of the EoS, 


r(r!) 


n dP 


(48) 
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Fig. 20. Neutron star masses (upper panels) and radii (lower panels) versus the 
central density for various EoSs without (left) and with (right) inclusion of the 
crust. Bold dots show maximum available masses. The band shows uncertainty 
range of the mass of PSR J0348+0432 (2.01 ± 0.04Mq). Horizontal dashes in the 
right panel show the DU thresholds. 

as a function of n varying from the star center to the surface (A), and the 
corresponding part of the star mass (C) accumulated at densities smaller than 
given n. Knowing the density prohle across the star, n{r), we present these 
quantities as functions of the radial coordinate in the right column in Fig. |2T] 
(B and D). Using the adiabatic index F as a measure of the local stiffness 
of the EoS we conclude that the cj-cut procedure stiffens strongly the EoS in 
comparison to the KVOR EoS in the narrow interval of densities above hq. 
Beyond this interval at higher density the KVORcut EoSs become even softer 
than KVOR EoS. This can be also observed on the accumulated mass plot. 
The accumulated mass M{n) for the KVORcut EoS remains smaller than that 
for the KVOR EoS at higher n inside the star and starts to exceed the mass 
for the KVOR EoS only in the outer part of the star. For the MKVOR EoS 
interval of densities, where its stiffness exceeds that of the KVOR EoS, is much 
broader than for any KVORcut EoS. As the result MKVOR EoS supports the 
neutron star with the larger maximum mass. 

In Fig. [22](left panel) we demonstrate the mass-radius relation for our EoSs 
in comparison with available experimental constraints. In contrast to the 
mass determination, there are no high-accuracy radius measurements. The 
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Fig. 21. Local stiffness (adiabatic index) of the EoS (A, B) and mass (C, D) distri¬ 
butions inside neutron stars with M = Mmax as functions of the density n (A, C) 
and the radial coordinate (B, D) for our EoSs. Crust is included. 



Eig. 22. Left panel: Mass-radius relation for our EoSs together with constraints from 
thermal radiation of the isolated neutron star RX J1856 |135] and from QPOs in 
the LMXBs 4U 0614-1-09 [133j . The band shows the uncertainty range of masses for 
PSR J0348-I-0432 [H]. Dashes indicate the DU thresholds. Right panel: Mass-radius 
relations for our EoSs together with other constrains on the star radius. Hatched 
regions show Icr (light), 2cr (hatched) and 3 cj (double-hatched) confidence contours 
from the mass-radius analysis [43] of the millisecond pulsar PSR J0437-4715. Thin 
lines show the 2a contour from BPA of Ref. [l6|. The elliptic region shows the 
estimated M-R constrained for RX J1856, cf. [136] . The rectangular box shows es¬ 
timations of compactness of the isolated neutron stars via X-ray spin phase-resolved 
spectroscopy from [44] . 
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highest observed quasi-periodic oscillation (QPO) frequency (1330 Hz, for 
4U 0614+091; see Ref. |133] ) places a constraint, see double-hatched area 
in the left panel. The nearby isolated neutron star RX J1856.5-3754 (shortly: 
RX J1856) shows a purely thermal spectrum in X-rays and in optical-UV. This 
allowed to determine the distance to the object |134j and its photospheric ra¬ 
dius Roo ~ 16.8 km |135] . Using the relation between the photospheric radius 
Roo and the true stellar radius R, Roo = R(1 — 2GM/R)~^/‘^ we can express 
the neutron star mass as a function of the radius as 

M = (2G')-i (^R - R^/Rl}j . (49) 

Assuming the estimate Roo ~ 16.8 km, which takes into account uncertainties 
in the distance and radius determinations, the hatched region in Fig. [22] Heft 
panel) is obtained as a constraint on the mass-radius relation |135j . There are 
three ways to interpret it: (i) RX J1856 is a typical single star with the mass of 
~ 1.4Mq, then its radius R must be larger than 14 km. Only our KVORcut02 
model satishes this constraint; (ii) RX J1856 has a typical radius of R =12- 
13 km, then its mass should be above (1.8-2) Mq, requiring EoSs like our 
KVORcut03 and MKVOR EoSs; (iii) RX J1856 is an exotic object with a small 
mass M ~ O. 2 M 0 , which would be possible for all EoSs considered here. Such 
small-mass neutron stars cannot be produced within standard models of the 
neutron star formation. Nevertheless, Ref. [138j suggested that fragmentation 
of rapidly rotating proto-neutron stars could lead to the formation of very light 
neutron stars. In a more detail these constraints are discussed in Ref. [31] . 

In Ref. [l3] the neutron star radius is constrained to be > 11.1 km at 3a 
conhdence, assuming M = 1.76Mq for all combinations of the other param¬ 
eters, see Fig. [22] (right panel). Thin lines in Fig. [22] (right panel) show the 
2a contour from Bayesian probability analyses (BPA) of Ref. jlB]. The el¬ 
liptic region ranging from 11 km to 12.7 km shows the estimated M-R con¬ 
straint for RX J1856, cf. [136] . We see that for all models, except KVORcut02 
model, the radii of the stars he in a narrow interval 12.5 ± 0.5km for masses 
O. 5 M 0 ^ M < 2Mq. This observation puts in doubt the hope expressed in 
some works that a simultaneous measurement of M and R may allow to £x 
the EoS with an appropriate certainty. The stars computed with the KVOR- 
cut02 model yield larger radii. The reason for this is that for the KVORcut02 
model the mass and the pressure are higher in peripheral star regions than for 
other EoSs that we consider, see Fig. [2TJ Note that with the help of X-ray spin 
phase-resolved spectroscopy. Ref. [H] found distinctly larger radii: i? > 14 km, 
as is shown by the rectangle. Among the hadron EoSs only the stiffest EoS 
for the KVORcut02 model may allow to fulhll this constraint. Perhaps, such a 
large radii occur for the so-called hybrid quark-hadron stars [139] . The reason 
why hybrid stars yield large radii is similar to that we have explained compar¬ 
ing our KVORcut02 EoS with other models: in hybrid star configurations one 
matches the hadron EoS in the shell with a softer quark one in the interior. 
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Thus, in order to satisfy the constraint M > 2Mq with hybrid stars one needs 
to exploit the stiffest hadron EoSs in the star shell. 

Because of mentioned ambiguities of evaluations of the neutron star radii we 
will consider as a relevant the constraint of Ref. |13]: R > 11.1 km at 3a 
conhdence for the star with M = I.IQMq. This constraint is well satished in 
all our models. 


6 Strangeness content. Hyperons, 0 mesons, aH scaling. 


In this section we study results of the hyperonization phase transition on the 
EoS within our KVOR-based models. The ratios of the aH to aN coupling 
constants deduced from hyperon binding energies in ISM following Eq. 0191) 
are: 


KVOR, KVORcut: = 0.599 , = 0.282 , = 0.305 , 

MKVOR; x^A = 0.607, x<,s = 0.378 , x^h = 0.307. (50) 

The values x^-a = 0.60-0.61 agree well with the best fit value derived from 
medium-heavy hypernuclei, x^a = 0.62, see Ref. m- 

After the inclusion of hyperons we label the KVOR, KVORcut and MKVOR 
models as KVORH, KVORHcut and MKVORH, respectively. Since with¬ 
out hyperons KVOR and KVORcutOd models yield rather low maximum 
masses and KVORcut02 model does not appropriately fulfill the flow con¬ 
straint for n < 4no, we will focus our further consideration on the best choices: 
KVORHcut03 and MKVORH models. As in above consideration, all calcula¬ 
tions below are performed with the BPS crust EoS. 

As known, the inclusion of hyperons results usually in a substantial decrease of 
the maximum neutron star mass. The difference between neutron star masses 
with and without hyperons proves to be so large for reasonable hyperon frac¬ 
tions in the standard RMF approach that in order to solve the puzzle one 
needs to start with a very stiff EoS without hyperons. This hardly coincides 
with the predictions of the microscopic variational EoS [32] and the quantum 
Monte Carlo simulations m, and is not compatible with the constraint on 
stiffness of the EoS extracted from the analysis of nucleon and kaon flows [3811^ 
in heavy-ion collision. 

In Fig. [23] (left panel) we show the dependence of the scalar field /(u), see 
Eq. ([8|), and the baryon fractions on the density for KVORH model (solid 
lines) and KVORHcut03 model (dashed lines) in BEM. In the right panel 
we demonstrate neutron star masses versus central densities for these models. 
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Fig. 23. Left panel: The scalar field amplitude /(n), Eq. ([8|), and the baryon 
fractions as functions of the total baryon density for KVORH model (solid lines) 
and KVORHcutOS model (dashed lines) in BEM. Vertical thin solid and dashed 
lines indicate maximally possible values of the central density for KVORH and 
KVORHcutOS models, respectively. Right panel: Neutron star mass versus central 
baryon density for the KVORH model (solid line) and for the KVORHcutOS model 
(dashed line). 

We see that the A and begin to appear for n ~ (2.5-3)no. Each of these 
fractions reaches ~ 0.2 for n > (5-6)no in both models nnder consideration. 
The hyperons appear only close to the maximnm available central density. 

do not appear at all. For the KVORHcntOS model the hyperon fractions 
prove to be larger than for the KVORH model. 


The maximnm mass is I.OOMq for the KVORH model and 1.83 Mq for the 
KVORHcnt03 model. Thus, the reductions of the maximum mass of the star 
equal to 0.35Mq and 0.33Mq for the KVORH and KVORHcut03 model, re¬ 
spectively, prove to be less than those would be found in the usual RMF 
approach [88]IWj . Nevertheless, the KVORH and KVORHcut03 models with 
the hyperon coupling constants given by Eq. fl50ll produce too low maximum 
neutron star mass. Moreover, we hnd that hyperons cannot be incorporated 
in the MKVOR model with coupling constants given by Eq. fl50|) and the rjp 
scaling function from fl29l) . since the solution for f{n) does not exist in this 
case. 


To have a possibility to increase the maximum values of the neutron star 
masses we incorporate the c^-meson mean held and allow for a scaling of the 
aH coupling constants. We will use the very same scaling of the (/)-meson 
mass as for all other hadrons, i.e., = 1 — /, but with the unsealed coupling 

constants, X(t>b = 1, which gives the scaling function = (1 — /)^. 

First, let us use the (pH coupling constants given by Eq. flTSD and the aH cou¬ 
pling constants given by Eq. flT^ with p,aH{f) = 1- So generalized KVORHcut03 
and MKVORH models are labeled as the KVORHcutO30 and MKVORH0 
models, respectively. 
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Fig. 24. The same as in Fig. [23] but for models KVORH(/)cut03 (dashed lines) and 
KVORH(/)(Tcut03 (solid lines). 

Then, additionally we incorporate 7^ 1 scaling. In order to demonstrate 
a qnalitative effect of the scaling we take ^cr//(n) such that ^^^(no) = 1, 
cf. Eq. (IT^ . and it decreases reaching zero for the baryon density n > ricH, 
iauin) = 0 for n > HcH) where Uch are critical densities for hyperonization. 
This assumption means that with such a parameterization we will exploit 
vacuum masses of the hyperon H ioi n > Uch- Such a models are labeled 
as KVORHcutO30(T and MKVORH^cr models. Note that the KVOR model 
extended to the high temperature regime in Ref. [35] (the SHMC model) 
matches well the lattice data up to T ~ 250 MeV provided all the baryon-a 
coupling constants except the nucleon ones are artificially suppressed, that 
partially motivates our choice of ^^h- 

On the left panel in Fig. [2l|we demonstrate f{n) and the baryon fractions for 
KVORHcutO30 model (dashed lines) and for KVORHcutO30cT model (solid 
lines). The particle fractions show a qualitatively similar behavior to that in 
KVORHcut03 model, except that now the do not appear at all and A do 
not occur for KVORHcutOd^a model. In the right panel we demonstrate neu¬ 
tron star masses as functions of the central densities for the same models. We 
get Mmav [KVORHcut03d)] = 1.97Mq that marginally agrees with the obser¬ 
vational neutron star mass constraint and Minax[KVORHcut03(;/)(T] = 2.08Mq 
that fully fulfills the maximum neutron star mass constraint. 

On the left panel in Fig.[25]we show f{n) and the baryon concentrations for the 
MKVORHc/) and MKVORH^cr models. The hyperon concentrations demon¬ 
strate a behavior similar to that in the KVORHcutO30 and KVORHcutOd^a 
models, being however suppressed by a factor ~ 0.7. On the right panel we 
demonstrate neutron star masses as functions of the central densities for these 
models. Mmax[MKVORH(;/)] proves to be 2.22Mq and Minax[KVORHcut03(;/)(T] = 
2.29Mq, both values are well above the maximum neutron star mass con¬ 
straint. 

Maximum values of the masses of neutron stars and the total strangeness 
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Fig. 25. The same as in Fig. [23] but for models MKVORH(/) (dashed lines) and 
MKVORH(/)ct (solid lines). 

Table 8 

Maximum masses, M^ax) and strangeness concentrations, /s of the neutron star for 
various RMF models with hyperons. 



KVOR 


KVORcut03 

MKVOR 


Afmax [Af©] 

/s [%] 

Afmax [AIq] 

/s [%] 

Ahmax [Af©] 

/s [%] 

no hyperons 

2.01 

0 

2.17 

0 

2.33 

0 

H 

1.66 

3.4 

1.84 

2.4 

- 

- 

H0 

1.88 

3.5 

1.97 

2.9 

2.22 

2.3 

H0cr 

1.96 

0.92 

2.08 

0.98 

2.29 

0.62 


concentrations calculated as the ratio of the number of strange quarks in the 
star to the total number of quarks given by 3Nb, where Nb is the baryon 
number, Eq. (06]), 

1 f dvrdrr^ 

/s = 17 ^ / ! ^=(nA + ns- + nso + ns+ + 2nH- + 2nso), (51) 

^l-2GM{r)/r 

for the KVOR, KVORcutOS and MKVOR models and their H0 and H0(T 
extensions are summarized in Table [3 We see that with inclusion of the 0 
meson mean-held with the scaled mass (H0 extensions of the models) the 
strangeness concentration increases but the maximum mass of the star in¬ 
creases too. Switching off the gaH coupling constants, as in H0cr extensions, 
suppresses /s drastically that leads to a further increase of Mmax- 

Hyperon critical densities and corresponding neutron star masses for KVORH, 
KVORHcutOS models and their 0 and (pa extensions and for extensions of 
the MKVORH model are presented in Table [9l We see that before the proton 
fraction reaches the DU threshold for n ^ p + e + u reaction, the DU reactions 
on hyperons, A p+e + 9, p+e A + u, and —)■ A-|-e-l-h, A-|-e —)■ S“-|-z/, 
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Table 9 

Hyperon critical densities in units of uq and corresponding neutron star masses in 
units of Mq for KVORH, KVORHcut03 and MKVORH models. Dashes mean that 
the hyperon species do not appear in the neutron star up to the maximal density. 



^cA 


ncE- 

[I] 

1 

O 

[1] 

O 

Mg' 

KVORH 

2.81 

1.37 

3.13 

1.48 

6.27 

1.66 

KVORHcut03 

2.59 

1.51 

2.89 

1.65 

5.10 

1.84 

KVORHcut03(() 

2.59 

1.51 

2.98 

1.69 

- 

- 

KVORHcut03(()(T 

- 

- 

3.42 

1.91 

- 

- 

MKVORH(/) 

2.63 

1.43 

2.93 

1.65 

- 

- 

MKVORH(/)a 

- 

- 

3.61 

2.07 

- 

- 



R [km] R [km] 


Fig. 26. Mass-radius relation for KVORHcutOSc/) (solid line), KVORHcut03i;A(T 
(dashed line), MKVORHc/) (long dash-dotted line) and MKVORHc/xr (dotted line) 
EoSs. Notations of the regions are the same as in Fig. \22[ 

may occur. As we see from the Table [9] in all considered cases except KVORH 
and MKVORH0, > l.SM®. Thus the KVORH and MKVORH0 models 
do not fulfill the “strong” DU constraint (Mdu > I-SMq), although all models 
satisfy the “weak” DU constraint (Mdu > I.SSMq). After one includes in the 
MKVORH0 model the < 1 scaling the “strong” DU constraint can be 
fulhlled, as we demonstrated on example of the MKVORH^cr model. Note 
that for KVORHcutO30(T and MKVORH^cr models reactions —)• A + e + u, 
A -|- e —?• -|- u are not allowed for n > Ucs-, since there are no A’s, and 

reactions n + e + v, n + e ^ 'E~ + v are forbidden due to the change 

of the strangeness by two units. Thus, in these cases values rics-, M^ mean 
the critical density and the mass for the appearance of rather than the 
corresponding DU threshold values. 

Mass-radius relations for KVORHcutO30, KVORHcutO30(T, MKVORH0 and 
MKVORH0CT EoSs are shown in Fig. [26](left and right). We see that all pre¬ 
sented models satisfy the maximum neutron star mass constraint. A general 


46 

































Fig. 27. The squared speed of sound as a function of the baryon density for various 
RMF models in the BEM with hyperons. 


behavior of the curves is similar to that shown in Fig. |22l where the hyperons 
and 0 meson are not incorporated, and = 1 - 


In the recent analysis [137] based on several theoretical EoSs of dense mat¬ 
ter the conclusion was drawn that neutron stars with masses I.O-I. 6 M 0 
with necessity have radius > 13 km. Our cut-extensions of the KVOR model 
with hyperons - KVORHcutO3(0cr) - do agree with this conclusion rendering 
radii ~ 13 km. However, the MKVORH(0cr) models have smaller radii ~ 12- 
12.5 km. The EoSs used in [137] are chosen very stiff in a pre-hyperon phase, 
i.e. at densities < 2-3no, by necessity, since they have to support a sufficiently 
heavy neutron star even after the softening of the EoS in the hyperon phase. 
The stiffening of the EoS at small densities leads to an increase of the star 
radius as we can see at examples of the KVORcut02 and KVORcut03 models. 
The latter two models and the EoSs used in [137] go above the constraint [3H] 
from the particle flow in HICs at least for n < 2 - 3 ? 7 ,o. Differently, the MKVOR 
model is made softer in this density range to pass the constraint [38] and is 
stiffened at higher densities to reach a sufficiently high maximum mass of the 
neutron star. This makes the star radius calculated in this model larger. 


Fig. 127] presents nf (n) as a function of the baryon density for BEM. The sharp 
bends appear at n = UcH) he. at the critical points of the third-order phase 
transitions associated with accumulation of the strangeness. For n < minjncff} 
the curves coincide with the corresponding curves in Fig. [TTl 


Thus the “hyperon puzzle” can be resolved in the framework of the models 
with scaled hadron masses and coupling constants provided 0 meson is in¬ 
cluded and its held mass scales at the same rate as masses of the a, uj and p 
helds but the (pH coupling constants are not scaled. 
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7 Conclusion 


In this paper we studied the equation of state of the cold nuclear matter con¬ 
structed within the relativistic mean-held approximation with hadron held- 
mass terms and coupling constants dependent on the a mean held variable 
/ oc o'Xo-Ar(o'), where X(Tiv(o') is the scaling function of the aN coupling con¬ 
stant. As the simplest choice, all hadron held masses are supposed to change 
with / at the same pace. The basics of the model were formulated in Ref. [2T] 
and a working model MW(n.u., = 0.65) labeled in Ref. [21] as KVOR has 
been constructed. In Ref. [3T| it was shown that the KVOR model appropri¬ 
ately satishes the majority of experimental constraints known to that time. In 
Refs. [3ni35ll36ll37] the KVOR model was generalized to non-zero temperatures, 
the particle excitations were included and the model was successfully applied 
to the description of heavy-ion collisions. Then Refs. [inBi] demonstrated new 
measurements of massive neutron stars, yielding the novel constraint on the 
maximum neutron star mass M^ax > 1.97Mq. The KVOR model, yielding 
Afmax > 2 .OIM 0 , hts this constraint, although marginally. 

The problem appears when one tries to include strangeness content within 
the standard relativistic mean-field framework. Owing to the hyperonization 
phase transition, the equations of state are strongly softened. Thereby it was 
concluded that only equations of state, being extremely stiff without hyperons, 
are able to explain observed most massive neutron stars. Inclusion of the 
interaction with 0-meson mean field, using artificially suppressed coupling 
constants, and other modihcations, as introduction of an additional repulsive 
terms into the Lagrangian, although soften the problem but do not fully solve 
it, see, e.g.. Ref. ra .Thus the aim of the given work was to show how one can 
construct an appropriate set of equations of state satisfying presently known 
experimental constraints within the relativistic mean-held models including 
hyperons and with hadron masses and coupling constants dependent of the 
scalar held. 

First, we considered models with an artihcially suppressed strangeness content. 
More specihcally we focused on two types of modihcations of the previously 
studied KVOR model. One type of models (labeled KVORcut) demonstrates 
that the equation of state stihens, provided f{n) < 1, if its monotonous in¬ 
crease for n ~ ritr (ritr > undergoes a sharp transition to a constant value. 
The smaller a transition value of / is chosen, at which the f{n) growth is 
quenched abruptly, the stiher becomes the resulting equation of state. On 
the example of the original Walecka model and a non-linear Walecka model 
we demonstrated that suggested cn-cut scheme is applicable to all relativistic 
mean-held equations of state, being a simple and efficient procedure that al¬ 
lows to stihen a given equation of state for n > Utr- In case of the KVOR-based 
models, this transition value of / is the smallest for KVORcut02 model and 
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increases for KVORcutOS and then for KVORcnt04 models. The sharp change 
of f{n) near Utr is provided by the change of the ratio of the scalings of the 
hadron effective masses to the Nuj meson conpling constant. 

The KVOR model exploits the effective nucleon mass m*{no) = 0.805m at 
at nuclear saturation and incompressibility K = 275 MeV. The new model 
(MKVOR) using the effective nucleon mass at nuclear saturation m^(no) = 
0.73mAT and incompressibility K = 240 MeV and newly adjusted scaling func¬ 
tions, gives the stiffer equation of state. In the MKVOR model the sharp 
change of f{n) near ritr in the beta-equilibrium neutron star matter is pro¬ 
vided by the change of the ratio of the scalings of the hadron effective masses 
to the Np meson coupling constant. 

We performed a comparison with results based on the microscopic equations of 
state such as the Urbana-Argonne group A18 -t- -|- UIX* one [32], the quan¬ 

tum Monte Carlo equation of state [96], the Skyrme parameterizations [65] . 
and the chiral symmetry based calculation ng. Also we calculated the Landau 
parameters /o, /i in the isospin symmetric matter and purely neutron matter 
and the parameters /g, f[ in the isospin symmetric matter as functions of 
the baryon density within our models and compared the results with those 
obtained with the Skyrme parameterizations. 

Then we confronted our results to the various constraints extracted from the 
numerous experimental data. More specihcally we compared the results with 
the nucleon optical potential, the heavy-ion collision nucleon and flow con¬ 
straints, the direct Urea constraint, the gravitational - baryon mass constraint, 
the maximum neutron star mass and mass-radius constraints and with some 
other constraints, as the symmetry energy constraint extracted from analysis 
of the isobaric analog states, the predictions of the neutron-proton elliptic flow 
differences confronted to the FOPI-LAND data, the giant monopole resonance 
data, etc. 

The most challenging is to fulhll the nucleon and kaon-flow constraints and, 
simultaneously, produce a high maximum mass of a neutron star. For that 
we introduced the scaling functions such that our equation of state is rather 
soft for n < 4no in the isospin symmetric matter but is sufficiently stiff in the 
/3-equilibrium matter. This behaviour is achieved by the proper selection of 
the scaling functions and Pp). 

As an interesting finding we indicate that the gravitational - baryon mass, 
Mq-Mb, constraint is better satisfied for models with the lower proton frac¬ 
tion in the density interval hq < n < 2.5no. The KVORcut02 model yields the 
largest proton fraction in the given density interval and does not fulfill the 
constraint. The proton fractions for KVOR, KVORcut04 and KVORcutOS 
models almost do not differ for n < 2.5no and their Mq-Mb lines almost coin- 
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cide in the considered Mq-Mb interval. Note that although these models have 
almost the same value of the density derivative of the symmetry energy per 
particle (h/3) as KVORcut02 model, the proton fractions of KVOR, KVOR- 
cut03, KVORcutOd models are essentially different from that for KVORcut02 
one. The model MKVOR is the best to satisfy the Mq-Mb constraint. In 
this model the proton fraction is the lowest among considered models in the 
density interval rio < n < 2.5nQ. Also we demonstrated that the fulhllment of 
the Mq-Mb constraint partially correlates with the value of L\ the smaller L 
is the better the constraint is satished. For L ~ 41 MeV even the “strong” 
constraint of Ref. |12] is fulhlled. 

Another result is that the KVORcut02 model yields a larger radius of the neu¬ 
tron star of the given mass than all other models considered since the energy 
density and the pressure in this model are higher at large radial coordinates, 
cf. Fig. i.e. for smaller densities n, than for other models. Since there 
are, although indirect, indications [43M4] on large radii of some neutron stars 
(R > 14 km), it might indicate that the equation of state that may describe 
these data should be very stiff for densities n < 4no and it can be softer for 
higher densities, e.g. such a softening might be associated with one of possible 
phase transitions at these densities, like the hadron-quark phase transition. 
However we should stress that with an additional stiffening of the equation of 
state for n < 4no the heavy-ion collision flow constraint is hardly satished. A 
weaker constraint of [43] that R > 11.1 km at 3cr conhdence for the star with 
M = 1.76M0 is well satished with all our equations of state. 

As the next step we included hyperons. As in all known schemes, accumulation 
of hyperon Fermi seas diminishes the maximum neutron star mass. However 
within our scheme this decrease of the mass is more moderate than in the 
standard relativistic mean-held approach. 

Then we incorporated the 0 meson mean held into consideration. We have 
demonstrated that with the inclusion of the 0 meson the maximum neutron 
star mass increases and the corresponding KVORHcutO30 and MKVORH0 
models satisfy the maximum neutron star mass constraint M > 1.97Mq, pro¬ 
vided 0 meson-baryon coupling constants are unsealed, but the 0 meson mass 
scales at the same rate as all other hadron held masses, Minax[KVORHcutO30] = 
1 .97M0 and Minax[MKVORH0] = 2.22Mq. Thus, the inclusion of the 0 me¬ 
son into consideration, such that the mass term scales at the same rate as for 
other hadron helds and the 0i7 coupling constants follow the SU(6) symmetry 
relations, allows to solve the part of the “hyperon puzzle” associated with a de¬ 
crease of the neutron star masses in the presence of hyperons. Then we demon¬ 
strated the ehect on the neutron star masses of the scaling of the aH coupling, 
taking it such that < no)) = 1 and iaH^fin)) decreases with 

n for n > hq. As the limiting case, we exploited the choice when the hy¬ 
peron masses reach their vacuum values for n > min{nciy}. In the correspond- 
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ing models named KVORHcutO30cT and MKVORH^cr, the maximum neutron 
star mass still increases, and we found Mmax[KVORHcutO30(T] = 2.07Mq and 
M^ax[MKVORH0a] = 2.29Mq. 

The low value of the critical density for the occurrence of A hyperons in the 
MKVORH^ model, Uca = 2.62no, may cause a problem with description of 
the neutron star cooling due to the opening of the efficient A —)■ p + e + z/, 
p + e —)■ A + z/ direct Urea reactions on the hyperon for M > = 1.43 Mq. 

This part of the “hyperon puzzle” is avoided within the extension of the model, 
if one allows for the scaling of the aH coupling. In the MKVORH^cr model 
A hyperons and hyperons do not appear. Thus, we are able to resolve the 
DU part of the “hyperon puzzle”. 

Other constraints considered in the given paper are satished equally good 
in the models with hyperons KVORHcutO30 and MKVORH(/), and KVORH- 
cutO30cr and MKVORH0cr, as it were in the models KVORcut03 and MKVOR 
in the absence of hyperons. 

In the given work we exploited traditional SU(6) based hyperon couplings, 
whereas we could use SU(3) ones, that would allow by varying of an extra 
parameter to further increase the maximum neutron star mass. We did not 
include the A isobars into consideration. Reference |140] raised a question 
about the A puzzle. They found that for 40 MeV < L < 62 MeV, A isobars 
in beta equilibrium matter may appear at a density of the order of 2—3 times 
nuclear matter saturation density, i.e. in the same range as for the appearance 
of hyperons. We believe that the A puzzle might be solved similarly to that 
for the hyperons. In a more detail this question will be considered elsewhere. 

A summary of equations of state considered in the given work in compari¬ 
son with the experimental constraints is presented in Table [101 We see that 
KVORcut03 and MKVOR models without hyperons and KVORHcutOO^a, 
MKVORH^a models with hyperons and </> meson pass the suggested tests 
rather appropriately. 


Acknowledgement 


This work was supported by the Ministry of Education and Science of the 
Russian Federation (Basic part), by the Slovak Grants No. VEGA-1/0469/15, 
and by “NewGompStar”, GOST Action MP1304. Gomputing was partially 
performed in the High Performance Gomputing Genter of the Matej Bel Uni¬ 
versity using the HPG infrastructure acquired in Project ITMS 26230120002 
and 26210120002 (Slovak infrastructure for high-performance computing) sup¬ 
ported by the Research & Development Operational Programme funded by the 


51 




ERDF. 


52 



Model 

^sym(^) 

IAS 

Fig. dSj 

HIC flow 

Fig. [T6| 

K~^ prod 

Fig. dSj 

DU 

w./str. 

Tables I6I9I 

Mg vs. Mb 

1%Mq 0%Mq 

Fig. dTj 

max. mass 

> 1.97Mo 

Fiffs.[22lMl 

star radius 

BPA J0437 

Fitrs. [22|26l 

KVOR 

o 

-b 

-b 

+/+ 

o 

- 

-b 

o 

+ 

KVORcut04 

o 

+ 

+ 

+/+ 

o 

- 

-b 

o 

+ 

KVORcut03 

o 

+ 

+ 

+/+ 

o 

- 

+ 

— 

+ 

KVORcut02 

o 

- 

- 

+/+ 

- 

- 

+ 

- 

+ 

MKVOR 

+ 

+ 

+ 

+/+ 

+ 

o 

+ 

+ 

+ 

KVORHcut03(/) 

o 

+ 

+ 

+/+ 

o 

— 

o 

o 

+ 

KVORHcut03i;A(T 

o 

+ 

+ 

+/+ 

o 

- 

+ 

o 

+ 

MKVORH(/> 

+ 

+ 

+ 

+/- 

+ 

o 

+ 

+ 

+ 

MKVORH0C7 

+ 

+ 

+ 

+/+ 

+ 

o 

+ 

+ 

+ 


Table 10 ^ 

Performance of various models proposed in text with respect to constraints on the nuclear symmetry energy (Tsym(^)) extracted from 
the analysis of the isobaric analog states (IAS) [67|; constraints on the pressure of the nuclear matter extracted from the analysis of a 
particle flow in heavy-ion collisions (HIC flow) [38] and a K~^ yield in heavy-ion collisions {K~^ prod.) |126j : weak and strong (w./str.) 
constraints on the thresholds of direct Urea (DU) reactions; constraints from the relation between gravitation and baryon masses {Mq 
vs. Mb) for the pulsar J0737-3039(B) [42| with account for the zero-mass loss during the proto-star collapse {0%Mq) and the loss of 1% 
of the Sun mass {1%Mq); constraint on the maximum mass of the neutron star (max. mass) [41] : constraint on the neutron star radius 
(star radius) extracted from the Bayesian probability analysis (BPA) of several neutron stars in [46] and from the analysis of the pulsar 
J0437-4715 (J0437) in |13|. Fulfillment (violation) of a constraint is marked by -|-(—) and a marginal result is marked by o. 














Fig. A.l. Pressure in the KVOR model matched with the BPS crust as a function 
of the number density (left panel) and as a function of the energy density (right 
panel). Solid line represents the resulting curve, dashed line is for the pure KVOR 
EoS, and long dash-dotted line is for the BPS EoS. 

A Matching the RMF EoS with the crust EoS 


The inclusion of the crust for all our models is performed as follows. There 
is always an intersection point between pressures as functions of the density 
in our hadron-lepton models for BEM and the BPS EoS frequently used to 
describe the neutron star crust [HI]. For KVOR and for all KVORcut models 
the intersection density is 0.47no, whereas for MKVOR model it is 0.65no. We 
have chosen a rather narrow interpolation interval for P{E), and above and 
below it we use exactly our EoS and the BPS P{E), respectively. The func¬ 
tion P{E) is what we need for integration of the Tolman-Oppenheimer-Volkoff 
equation to obtain the star mass and radius, which are of our main interest 
(especially M). The interpolation interval in E, that we have chosen, corre¬ 
sponds to the density interval 0.45no < n < 0.7no- Interpolating EoS in this 
density interval, we use a cubic spline. To get thermodynamically consistent 
P{n) and E{n) dependencies we reconstruct the baryon number density n as 
a function of E using the equation 



ni 



n 


dE 

P{E) + E' 


(A.l) 


which we integrate from a density rii > 0.7 uq downwards to lower densities. 
We chose rii above but not far from 0.7rio, so that at n ~ rii all KVORcut 
models are the same. For n > 0.7?7,o we preserve our model EoSs exactly. Here 
we follow the logic similar to that in Ref. |142| . As the authors indicate there, 
such a procedure in not unique. However, for the present calculations, where 
the main ingredient for us is the function P{E), it seems us optimal. The 
so-reconstructed function P{n) coincides with P{n) from our model (KVOR, 
KVORcut or MKVOR) for n > 0.7no and agrees well (on the fraction-of-a- 
percent level) with P{n) of the original BPS EoS for n < O.dhrio- Thus, the 
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E{n) dependence is thermodynamically consistent with P{n) and P{E) at all 
densities. The resulting functions P{n) and P{E) are shown in Fig. lA.ll for 
the KVOR model. The EoSs for the KVORcut02, KVORcutOS, KVORcut04 
models are exactly the same as for the KVOR model in the given density 
region. 


The region of densities, where we exploit interpolating EoS, is actually the 
region of a pasta phase, see |143] and references therein. The presence of the 
pasta affects EoS only slightly. Therefore, simplifying the consideration, we 
disregard this complication in this work. Note that the BPS crust EoS was 
also used in [31], where it was joined with the various EoSs describing the 
interior region. We have also checked that a narrowing of the spline-interval 
almost does not reflect on such observables, as the star mass and radius. 


B Derivation of Landau parameters 


Derivations of the Landau parameters /o and fi are performed following the 
lines of Ref. m- The parameter /q can be found with the help of our energy 
density functional flT^ exploiting Eqs. - fl3^ . In order to hnd the Landau 
parameter /i we need to keep the vector u and terms in the energy density. 


Here we present resulting expressions for the Fermi liquid parameters /o and 
/i, which we calculated for our generalized RMF models with scaled hadron 
masses and couplings for ISM and PNM. 

In case of the ISM (7 = 4 ): 


/o = 








Cln dr]^ m 
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^df 


X 


/i = 


:Y’fo(/)) 
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mjsf 


df 


+ 




4 U"{f) + 
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ClClf 1 

2m% \r]M) 


pI 


Pf + 


l + MnJ) 
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(B.l) 

+ Ki{nJ) + K2{n, f)J , 

(B.2) 
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^ J w^dp lp^ + m*j^ 

^nVM) { 27r2 (p2+|7^^)3/2 

Pl 

3m%p^{f)Tl^ ^p2 ^ ^*2 ’ 


(B.3) 


/(n) is solution of Eq. flT^ . Prime and double prime denote first and second 
derivatives over /, respectively. 


After the replacement 

Cl Cl Cl 

vM) ^ vM) ^ %(/) ’ 

expressions (140]) . (14T]) and fIB.ip . flB.2p become valid for the PNM (7 = 2). 
For the isospin-dependent Landau parameters in the ISM we get: 


foin) 

f[{n) 


Cl 

4m^7p(/) ’ 

Cl Pi 

^'mliPp{f) p], + m 


*2 

N 




(B.4) 

(B.5) 


KiinJ) 


Cl J ^p^dp lp^ + m*jl 

4m^7p(/) i 27r2 (p2 _|_ ^*^'j3/2 

7 Cl _ Pf 

8 3m^7p(/)7r2 ^^2 ^ ^*2 ' 


(B.6) 


Another way is to find the Landau parameters /o, /i and /q using expressions 
m. dUD and fl42l) . We have checked that with expressions (E^, flB.2p for 
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/o, /i, and dEl) for Eqs. 


fHT]) and (H2|) are indeed fnlfilled. 
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